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We consider the quantum dynamics of a magnetic domain wall at low temperatures, where dissi- 
pative couplings to magnons and electrons are very small. The wall motion is then determined by 
its coupling to phonons and nuclear spins, and any pinning potentials. 
I In the absence of nuclear spins there is a dominant superOhmic 1-phonon coupling to the wall 

velocity, plus a strongly T-dependent Ohmic coupling to pairs of phonons. There is also a T- 
independent Ohmic coupling between single phonons and the wall chirality, which suppresses "chi- 
rality tunneling. We calculate the effect of these couplings on the T-dependent tunneling rate of a 
^ • wall out of a pinning potential. 

'^'n ' Nuclear spins have a very strong and hitherto unsuspected influence on domain wall dynamics, 

1^-^ , coming from a hyperfine-mediated coupling to the domain wall position. For fcsT ^ ujq this 

' coupling yields a spatially random potential, fluctuating at a rate governed by the nuclear T2. 

When fcsT <^ luq, the hyperfine potential fluctuates around a linear binding potential. 

The wall dynamics is influenced by the fluctuations of this potential, ie., by the nuclear spin 
dynamics. Wall tunneling can occur when fluctuations open an occasional "tunneling window". 
This changes the crossover to tunneling and also causes a slow "wandering", in time, of the energy 
levels associated with domain wall motion inside the pinning potential. This effect is fairly weak in 
Ni- and Fe- based magnets, and we give an approximate treatment of its effect on the tunneling 
dynamics, as well as a discussion of the relationship to recent domain wall tunneling experiments. 
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S ; I. INTRODUCTION 

Ch The quantum dynamics of mesoscopic or even macroscopic magnetic solitons has been the subject of considerable 

^ interest in the last few years. There are several reasons for this. One is the rapid developraent of the new field of 
. . , nanomagnetismlJ, in which very small magnetic wires, films, and magnetic macromoleculesH (as well as molecular 
^ ■ chains) can be prepared. In many cases one expects that solitons such as domain walls, or perhaps lower-dimensional 
vortices and skyrmions, will play an important role in the quantum dynamics of nanomagnets at low temperatures in 
systems as diverse as magnetic perovskites, the cuprate superconductors, or nanomagnetic wires, films, and junctions. 
' The magnetic dynamics also strongly influence the charge transport in some systems. 

A second reason for this interest stems from the continued fascination with macroscopic quantum phenomena. 
Previously associated almost exclusively with superfluids and superconductoes, it has become clear in recent years 
that magnetic systems may also exhibit macroscopic quantum phenomenaU'E , -perhaps the most spectacular being 
the possible quantum tunneling of macroscopic domain walls-jB. ferromagnetso'El. There now exists some experimen- 
tal supnort_for tunneling of single domain walls in Ni wirea3 as well as many earlier experiments on multi-wall 
systemsE2rt3. The tunneling dynamics of large fiexible domain walls poses many interesting theoretical problems, 
particularly concerning the role played by the various environmental degrees of freedom, which couple to the wall in 
various ways. 

This is the main subject of this paper, which deals with the hitherto neglected effects of phorwHS and nuclear 
spins on magnetic domain wall tunneling (a limited discussion of phonon effects was given previouslyQ'El, but was very 
incomplete). We find that the effect of the nuclear spins in particular can be quite spectacular, and can completely 
contradict the picture of quantum domain wall dynamics that has evolved over the last few years. We derive a new 
scenario, in which magnetic domain wall dynamics is influenced at low temperatures by the effective coupling of the 
wall to a random (in space) "hyperfine field", which also fluctuates in time. Even in Ni wires (where these effects 
are far weaker than in other magnets), this fluctuating hyperfine field will be shown to have some influence on wall 
tunneling. The phonons are found to have a much weaker effect, but one nevertheless important to evaluate when one 
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considers the finite-T tunneling rate (particularly in isotopically purified systems, where phonons entirely determine 
the low-T dissipative dynamics, if the system is pure and insulating). We shall also see that phonons couple very 
strongly to the domain wall chirality, and thereby suppress chirality tunneling. 

In the remainder of this section we give a brief explanation of previous work in this field, and note how the present 
paper fits into this. 



A. Macroscopic Quantum Phenomena in Magnets 

The history of macroscopic quantum phenomena goes back to the 1930's (Meissner effect, superfluid fountain 
effect, capillary superflow, etc), and has until recently been associated almost exclusively with superconductors and 
superfluids. In recent years considerable attention has focused on Caldeira and Leggett.'aj)roposal that SQUIDs 
and related systems might show macroscopiG-|flux tunnelingtJ. Subsequent experimentsEjllj confirmed their work 
quantitatively, and this is still an active ficldlla. 

It is obviously important to find other systems in which quantum macroscopic phenomena occur. It should be 
stressed that this is not just a question of energetics but also involves the "environmental decoherence" which al- 
most always destroys quantum phase correlations at all but the microscopic level. Consider, eg., the possibility of 
"macroscopic quantum tunneling" (MQT) just mentioned. Many authors have suggested MQT after having noted 
the existence of some small tunneling barrier afitiag on some large coUertive mode - there i.s a long list of candidates 
for such MQT. The list includes ^He superfluidEj'Elland "^He superfluidcHl, superconductorst§E3 ferroelectricsEj, charge 



density waves (CDW) tunneling iijuJ-dimensional systcmsO, MQT and tunneling nucleatioii-of dislocations, vacanc 
and impuritiea-4n quantum solidsE3, quantum diffusion of such objects in quantum solidsEH, as-well as in norma 
and superfluidE3 ^He. To this-list we may also add more general studies of quantum nucleationE3 as well as of large 
scale tunneling in cosmologyEil This list in certainly not complete. Finally there are the various kinds of macroscopic 
quantum phenomena, including MQT, which have been suggested for magnetic systems, to which we come below. 

Despite the length of this list, in few cases have these suggestions been confirmed by experiments. There are 
two main reasons for this. First, the the surrounding "quantum environment" usually strongly suppresses MQT 
(aiuLcoHipletely suppresses "macroscopic quantum coherence" (MQT)). Consequently, as emphasised by Leggett et. 
aOEdLJ any discussion of MQT which ignores the environmental effects is usually of of academic interest only. 

A second reason is simply that most physical materials contain defects, impurities and other imperfections. The 
coupling of these to the collective coordinates involved in MQT is often much larger than the very small tunneling 
barriers. Since such external effects are usually unknown, this makes experimental tests very difficult. 

So far, the only case where really indisputablej-evidence can be given for the observati on ! O if MQT is in supercon- 
ductors (in SQUIDs and other Josephson devicestS). This work has been f*44^1y reviewedt§EJ. In magnetic systems, 
suggestions that tunneling phenomena might exist go back over 40 yearsLjLJ. These suggestions typically involved 
some 1-dimensional phcnomenological potential (a "pinning potential" for domain walls, an anisotropy potential for 
monodomain grains) to which a naive WKB analysis was applied. No microscopic theory was attempted, and the 
experimental evidence lacked any clear theoretical basis. 

However, the picture that has emerged in recent years is very different from this, for 2 main reasons. The first is the 
realisation that a standard WKB analysis is inappropriate to spin systems, in which a constraint typically exists - in 
small particles or macromolecules at low T the magnetisaticm density |]V[| is roughly constant (although not necessarily 
in domain walls, if multi-magnon excitations are importanttl) . This constraint leads to a kinetic energy which is linear 
in time derivatives (not quadratic^) J^hus a correct solution for the problem of tunneling of a spin was not given until 
1986, by van Hemmen and SutoEZrEa (and independently for a particular case by Enz and SchillingEJ) . Both of these 
analyses were semi-classical. A key role in nanomagnetic dynamics is also played by the .."Kramers/Haldane phase", 
which differs for integer and half-integer spins; this has no analogue in superconductortSo'cS. 

The second reason for the revised picture is that once microscopic theories were attempted of domain wall tunneling, 
or of tunneling in nanomagnets, it became obvious that yet again, much of the physics was in the coupling to 
the background quantum environment. Some of this physics is very different from that in superconductors. In 
particular, the coupling to nuclear spins is very strong, can play a major role in controlling the quantum dynamics 
of mesoscopic orpOEppscopic magnets, and cannot obviously be understood in terms of an "oscillator bath" model of 
the environmen10L3. Extensive theory has been done in an effort tOjimderstaiid the way in which the environment 
determines nanomagnetic dynamics, both in particles and moleculesLiLJuZI'HEj and in domain wall dynamics (see 
below). 

Since this theoretical work was done, strong experimental evidence has emerged for tunneling phenomena in meso- 
scopic magnetic systeros-at low temperature. This includes many experiments on multi-grain systema£j and on systems 
with multiple domainst^; evidence has been claimed both for magnetisation tunneling in monodomain systems and for 
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domain wall tunneling in puilti-domain systems. There has also been a claim for the observation of coherent tunneling 
in large ferritin maleculeS, which has been disputed by various authorg£j. A very interesting set of experiments by 
Giordano et. alflia reports strong evidence for tunneling of single domain walls in Ni wires, to which we will return 
below. Finally, a whole series of experimenis, on Mni2-acetate molecular crystals gives strong evidence for thermally 
activated resonant tunneling in this systemEJ. 

B. Domain Wall Tunneling 

The principal results to emerge from the first microscopic analyses of domain wall tunneling have been: 

(a) The effective pinning potential for a wall is almost independent of the nature of the pinning centre, provided the 
spatial extent of the pinning centre (usualljusome kind of defect) is much smaller than the wall width A; the shape of 
the defect is also not important in this easel. The pinning potential is of the form V{Q) = — Vosech^(Q/A), where Vq 
is positive. Subsequent analyses of the so-called "domain wall junction" rely on this result, which is important for any 
comparison with experiment, simply because the exact nature of the pinning centre is not known. The constant Vo 
can be determined experimentally in various ways (eg., from low frequency "wall rocking" absorption measurements, 
or, less reliably, form a determination of the wall escape characteristics at higher T). 

(b) For an insulating magnet, the principal contribution to wall dissipation at temperatures above the magnon gap 
A will come from multi-magnon emission by the walla. However if kT ^ A, the main magnon contribution to the 
dissipation will come from processes involving 2 bulk magnons and one wall, or "Winter" , magnon; this 3-magnon 
process gives an Ohmic contribution to the wall dissipation coefficient rj, going like 773 ~ {kT/ A) exp(— A/A:T). Higher 
magnon processes give contributions 77™ ^ (fcr/A)'"~^ exp(— A/ZcT), and are thus subdominant. 2-magnon processes 
give a superOhmic contribution to the wall dissipation which is negligible compared to the Ohmic contribution at low 
energy or low wall velocity. 

(c) When kT <C A, so that the magnona have negligible effects on the wall dynamics, an analysis has been given of 
2 other small contributions to dissipation^: a coupling between the wall and Winter magnons, mediated by defects; 
and also a quite negligible coupling, to photons. The influence of the defect-mediated coupling has recently been 
re-examined by Leduc and StampE3, with the conclusion that it can usually be neglected in practically relevant 
calculations. i_. 

(d) For a conducting magnetic system, a thorough analysis by Tatara and FukuyamaEJ showed that both the charge 
fluctuation (eddy current) and spin fluctuation contributions to the wall dissipation will be negligible unless the wall 
thickness is not much larger than the lattice separation between magnetic ions (even for SmCos, where A = 12A, the 
electronic contribution gave a correction < 1% to the tunneling exponent). 1-1 m 

(e) The_phonon contribution to the dissipation of wall motion was analysedH'B using the method of Wada and 
SchriefferEa to give a small T^ contribution to the wall diffusion constant. This analysis was incomplete, since it did 
not include all one-phonon contributions to the dissipation; moreover, it did not include higher non-linear couplings, 
such as 2-phonon couplings. Since the 1-phonon coupling is superOhmic, experience with magnon couplings suggests 
one should go to higher order to look for Ohmic contributions to the wall dissipation. This is one of the main tasks 
of the present paper. rnn 

(f) One can also consider the possibility of "chiralitft tunneling" of the domain wall chiralityEj'E2l, or even the 
coherent Bloch dynamics of walls in a periodic potentialEiJ. Until now no attempt has been made to see whether these 
possibilities would survive the coupling to the environment. This issue is also addressed herein; we find that they do 
not survive. 

Amongst the various experimental [Searches for macroscopic domain wall tunneling, perhaps the most dramatic 
results are those of Giordano et. alJsia. These results are in fairly good agreement with the theory, except that the 
experimental result for the quantum/classical crossover temp|erature (between tunneling and activated escape) is over 
an order of magnitude higher than the theoretical predictionO. 

In this paper we will show that yet another modification of our understanding of the low-T dynamics of domain walls 
(and other magnetic solitons) is necessary. This is because of the inevitable coupling of the wall to the background 
nuclear spins, via hyperfine interactions. We shall see that the longitudinal hyperfine coupling energy ^i-i^ksf.!^, 
summed over all nuclear spins 1^ and electronic spins s^^) in the system, can be very large. The transverse hyperfine 
coupling 'Ylik^k{s'^Ik + ^k-^k) causes irreversible transitions in the nuclear spin bath every time the nanomagnetic 
system changes its state. This changes the energetics of the tunneling, and also causes "topological decoherence" 
and dissipation. These names are imported from the study of environmental "spin bath" effects on nanomagnetic 
tunnelingc3u3, where similar physics applies. 

We will not try here to give a complete analysis of nuclear spins effects on domain wall dynamics - this would 
be a very lengthy task. We concentrate instead on the static longitudinal hyperfine coupling. The main conclusion. 
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developed in detail below, is that domain wall tunneling is rendered much more complex by the coupling to nuclear 
spins and to phonons. The longitudinal hyperfine field is random in space, and fluctuates in time at a rate governed 
principally by T2 relaxation. It can be very strong, and typically pins the domain wall. Tunneling can occur when 
temporal fluctuations in this potential open up "energy windows" for brief periods of time. The phonons also allow 
inelastic tunneling processes. Analysis of the various couplings to phonons shows that there are 3 important ones: a 
linear coupling to the domain wall velocity, a bilinear Ohmic coupling to the wall position, and an important linear 
Ohmic coupling to the wall chirality, which strongly suppresses any chirality tunneling, even in the absence of nuclear 
spins. As might have been expected, Bloch coherence is so fragile that it is destroyed by almost anything, and we 
believe it to be practically unobservable. 

The plan of this paper is as follows. In the next section we discuss the various domain walls we will be dealing with 
in various geometries, including magnetic wires, magnetic films, and bulk 3-dimcnsional magnets. Essential results 
for the wall dynamics and tunneling in the absence of environmental effects are given, as well as relevant experimental 
numbers. 

In section III we derive the effective low-energy interactions which couple a domain wall to both phonons and 
nuclear spins. The wall is coupled simultaneously to an "oscillator bath" of phonons, and a "spin bath" of nuclear 
spins, and a number of general results for such baths are recalled. 

In section IV an extensive analysis is given of phonon effects on domain wall dynamics, for walls in 1, 2, and 3 
dimensions, in isotopically purified magnets (no nuclear spins). The effective action for the wall is derived, with 
the phonons integrated out. These calculations are interesting theoretically, because of the novelty of some of the 
couplings. The dominant couplings to the wall velocity, the wall position, and the wall chirality are treated separately. 

In section V the problem of wall tunneling is analysed. We first consider the "displacement tunneling" of walls 
in isotopically purified magnets, where it is controlled by phonons. We then include the nuclear spins, and show 
that tunneling is then strongly influenced by the dynamic hyperfine potential acting on the wall. For displacement 
tunneling this changes the theoretical predictions for the tunneling (in particular the crossover temperature and the 
tunneling rate). As a corollary to these calculations we analyse the suppression of both chirality tunneling and of 
Bloch coherence. 

In the final section we make a few remarks on the relation to experiments on wall tunneling (particularly those of 
Giordano et al.). 

II. DOMAIN WALLS AND THEIR ENERGETICS 

In this section, we introduce the different kinds of domain wall we wish to study. Discussion of the couplings to 
the environment is reserved to later sections. We work exclusively with ferromagnetic systems, at low T, and we 
assume that the domain wall thickness A is sufficiently greater than the lattice separation ao between spins so that a 
continuum approximation for the magnetisation is valid. We will be discussing domain walls in non-trivial geometries, 
including films and wires, as well as 3-dimensional systems. 

In all cases, we assume that the magnetisation at position r in the crystal is defined as M(r) = g^B J2j 
where hb is the Bohr magneton fj,B = eTi/2meC, g is the Lande factor, and R-,- is the position of the j*'' lattice site, 
with associated electronic spin Sj . 

There is no easy way to deal with the magnetic structure in bulk materials and in thin films and wires in a 
unified manner, as the direction of the magnetisation within the walls is quite different in these two cases. Instead 
of considering different values for the magnetisation in the (x,y,z) directions, we define in each case (bulk and 
thin films) a new frame of reference (xi, X2, X3) such that the easy axis is along the xi direction, with the wall 
characterised by its position along the X3 axis. The magnetisation is written as M = Morii with the unit vector 
rfi = {'mx,'my,'rnz) = (cos sin 0, sin sin 6*, cos 0) so that (Vm)^ = (V6')^ -I- sin^ 6'(V0)^. We will assume that the 
magnitude Mq is constant, since all phenomena considered occur at temperatures far below the Curie temperature. 
We thus use A/q — sjgh/ag where 7^ — g^B/^ is the gyromagnetic factor, ao is the lattice spacing of the crystal and 
s is the value of the electronic-spin. We will assume a standard biaxial (easy axis/easy plane) Hamiltonian, given in 
continuum approximation byl£3 

W = i y dr[J(VM)2 - if||(Afi)2 + K^iM^f - ^(H^™ + H,,) • M] (2.1) 

where (VM)^ = (ViAfj)^. The easy axis is represented by the component Mi of the magnetisation, and the easy 
plane is perpendicular tn-M^. K^^ and K± are the anisotropy constants, He is an external magnetic field and Udm is 
the demagnetising fielcS, the internal field coming from all the magnetic moments. 
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The action of a [the magnetisation also contains a term with no classical analogue, the Berry phase factor0, 
Wess-Zumino termEj, given by 



l/T 



dr0(r,t)(l - cos6l(r,t)) 



(2.2) 



so that the total action is 5* = Swz + / drTi.. The Wess-Zumino term is purely imaginary. It essentially corresponds 
to the total solid angle traced by the spin configuration for a given trajectory; its derivation using path-iajtegral or 
coherent spin state descriptions of the system is now very well known, following the 1983 work of HaldandHj. 



1. Block Wall in Bulk Materials 



For our study of the tunneling of a Bloch wall, we take the Hamiltonian 
n = j dr[J{Vraf - K\\ml+ Ki_ml] 

= j dr[J{{\/ef + sin^ e[V<l)f) - K\\ cos^ e + Kx_ cos^ ^sin^ 



(2.3) 



(2.4) 



representing a ferromagnet with easy z axis and easy z — y plane. In 3-dimensions, J is measured in J/m, and the 
anisotropy constants in J/m'^. 

The domain wall corresponding to this Hamiltonian is perpendicular to the x axis, with the magnetisation rotating 
in the z — y plane. The wall centre is located at a position Q along the x axis. The new frame of reference is thus 
(xi,a;2,X3) = {z,y,x). This is represented in Fig.([|).|— , 

The components of the magnetisation are given byE3: 

=Ctanh(^i^^) 
.-.(l-||))sech(^) 

mi = C^sech 
^ 2co \ Xb J 

where the thickness of the wall \b — {J/K^ Y^^ represents the usual compromise between exchange and anisotropy. 
The surface energy of this type of wall is ao = 4(J/C|[)^/^. C = ±1 is the "topological charge" of the wall and 
X = ±1 is the "chirality". The topological charge corresponds to the direction along which the wall moves under the 
application of an external magnetic field in direction parallel to the easy axis, while the chirality refers to the sense of 
the rotation of the magnetisation inside the wall. The internal configuration of the wall is such as to simultaneously 
minimize the anisotropy energy (coming from both and Kj_), and the exchange energy. The anisotropy K± can 
originate from the material itself but also from the configuration of the magnetisation, through the demagnetisation 
energyES, ie., K± = K±,i + /ioMg/2 where K±^i is the anisotropy intrinsic to the material. A static Bloch wall only 
rotates in the easy plane. However, as soon as it moves it creates a demagnetising field which causes the spins to 
precess .Find .the appearance of a component of the magnetisation out of the plane, directly proportional to the wall 
velocityE^lEll. This picture is valid provided that Q{t) « cq, the Walker critical velocity, whereEj 



Co 
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(2.5) 



where typically cq ~ 10^ m/s. The precession of the spins also causes_the appearance of an inertial term, the Doring 
mass, given by the ratio of the wall energy with the limiting velocit5i£3: 



(2.6) 



where Sw is the surface of the wall. We assume K±^ ~ iiqM^/2 ^ and in this limit, cq — /io7gAMo/2 and the 
Doring mass reduces to — 2Sw/ fJ-olq^- 



5 



All of the above formulae for the Bloch wall are well known, and should be treated as phenomenological; they are 
given in terms of microscopic parameters, but should not be treated as the result of a microscopic derivation. In 
fact such a derivation would certainly give somewhat different results in the long-wavelength limit appropriate to the 
macroscopic dynamics of the wall, since integrating out the short wavelength fluctuations in M(r,t) will renormalise 
all parameters. We do not bother to do this for 2 reasons, viz., 

(a) It is obvious that the renormalising influence of the high-frequency fluctuations will be strongly suppressed by 
the long-range demagnetisation field 

(b) The only important parameter in any experimental tests of the theory will be the Doring mass, and this will have 
to measured anyway (by, eg., looking at the fluctuation frequency of the wall in a pinning potential); the measured 
mass will of course then already incorporate the renormalisations. 



Because of (a) we expect that the expressions in (2.5) and (2.6) to be fairly accurate anyway, and we will use them 
when necessary. 

If defects are present in the sample, they can pin the wall. We will assume that the radius Rd corresponding to the 
defect volume is much smaller than As, the domain wall width. The wall is then pinned by a potential of formQ 

ViQ) = ~Vosech\Q/\B) (2.7) 

with Vq proportional to the volume of the defectS. We further assume that there is a very small concentration of defects, 
so there is only 1 important pinning centre for the wall. This would correspond to an ideal experimental situation. 
We also assume that the wall is flat and that it remains flat during the tunneling process. This approximation is 
justifled by the energy associated with the curvature of the wall, to which there are two contributions. The flrst is the 
surface energy ctq of the wall. A curved wall has a larger surface, and thus a larger energy. Second, and much more 
important, a curvature in the xi — plane creates strong long-range demagnetisation flelds, which rapidly increase 
the wall energy. A complete treatment of these effects is quite complicatedE3, but for small pinning energies, the 
radius of curvatumjaf the wall is much larger than A, and it is easily shown that weak curvature has very little effect 
on wall tunnelingpEa. To consider the effect of dissipation, we can thus use a flat wall. 

The application of an external magnetic field He in the direction of the easy axis couples to the magnetisation to 
give a potential term linear in Q. Including the incrtial mass, we then write a "bare" Hamiltonian (ie., neglecting the 
environment) for the walfi 

Hw - ^M^Q^ - ViQ) - 2S^fiaMQHeQ (2.8) 

where we put the chirality C ~ 1 for brevity. 

The tunneling rate of the wall can now be evaluated by standard instanton techniques. For a field sufficiently close 
to the coercive field, the "pinning potential plus field potential" reduces to a quadratic plus cubic potential, whose 
barrier height is controlled by e = 1 — He/ He, He being the coercive field, given as 

while the energy barrier is 

V{e) ~ {h^gfioHe)Noe^/^ (2.10) 

where Nq is the number of spins in the wall. The exact value of the coercive field is of course very difficult to obtain 
theoretically and should ideally be obtained by a characterisation of the suateoi in the thermal phase. 

The tunneling rate Fq in a cubic plus quadratic potential is well knowroOllj; we simply recall here the results for 
the tunneling of the wall in the absence of dissipation: 



30B(e) 



1/2 

Qo e-^"'^"^/'' (2.11) 



with Bo{e) the usual WKB tunneling exponent; 



^Boie) = ^M^noQl (2.12) 
n 15 



flo the oscillation frequency of the wall in the potential; 
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(2.13) 



Qo is the escape point; 



(2.14) 



and A is again the width of the waU. In terms of the microscopic parameters, Bo{e) can be rewritten as 




54 5"^ A 
5 -fgh 



) 



1/2 



Mo 



C 



(2.15) 



where the second form comes from using the number of spins in the wah, Nq — \S^^/a^. 

Although these results were originally derived for the particular casef*)£ a short-ranged defect pinning potential, 
they obviously have more general applicability; as emphasized previouslyufl, for small e almost any pinning potential 
acting on the wall via the dipolar field or the wall surface energy will have the cubic-quadratic form near the coercive 
field, since its range will be ~ A. However we shall see later in this paper that the same is not true of the longitudinal 
hyperfine field, coming from the nuclear spins, which fluctuates over much shorter length scales. 

Let us at this point introduce 2 examples to which we will return at various points in this paper. We consider two 
particular systems, Yttrium Iron Garnet (YIG) and nickel. YIG is an insulator with a bcc cubic structure, a saturation 
magnetisation iiqMq = 0.24 T and with exchange and anisotropy energies J = 1 x 10^^^ J/ni and K^^ — 580 J/m . The 
width of the domain wall is A = 860A, with a mass per unit area 2 x 10~^ kg/m^. Nickel is a conductor, again with a 
cubic structure. The saturation magnetisation ^qMq = 0.6 T, with exchange and anisotropy J = 3 x 10~^^ J/m and 
ii'll = 4500 J/ni"^, giving a domain width A = 500A and a mass 6 x 10~^°kg/m^. As an example, we give the tunneling 
rate of a wall containing iVo — 10^ spins in nickel. Assuming e = 10^^ and Hc/Mq = 0.01, we get Bq/Ti ~ 20, and a 
frequency fig ~ 6 x 10^ sec~^ so that F ~ 200sec~^. 

The crossover temperature Tq between thermally activated relaxation and and quantum tunneling is roughly k^TQ ^ 
Q.q/2tt. For many potentials (including the quadratic/cubic potential) there is a fast crossover between the two modes 
of relaxation (there are however systems, such as Mni2Ac, where tunneling apparently takes place at intermediary 
levels, over a wide temperature range). For the Ni wall described above, Tq ~ 0.02 K; a crossover temperature 
To ~ O.IK would require the product /io-ffcC^^^ 0.6 T, ie., a coercive field similar to the magnetisation density. 



Bloch walls occur predominantly in bulk materials, and so are difficult to observe individually. Furthermore, they 
are fairly big, making the observation of tunneling difficult. If one tries to reduce the size of the sample to a wire or a 
platelet, then the Bloch wall becomes unstable and the magnetisation profile becomes quite complicated, most of the 
time being two-dimensional. One alternative is to go to thin films, of width 5 <^ \, and with an anisotropy axis in the 
plane of the film. In this case, the strong demagnetisation field forces the rotation of the magnetisation between, two 
domains to also take place in the plane of the film, in a head- on fashion. This is the model of the Neel wallElrEj. Its 
energy characteristics (ie., the Walker limiting velocity and the Doring mass), as well as the magnetisation profile are 
similar to that of a Bloch wall. We consider a film of length L ^ A, width 5 <C A in the xi — plane and with easy 
axis in the xi direction. The frame of reference of the wall is related to the Cartesian frame as (xi, a;2, 2:^3) = (y, z, x). 
This is shown in Fig. (^. 

The continuum Hamiltonian is 



2. Neel Walls m Thm Films 




(2.16) 



The profile of the magnetisation is thus 




(2.17) 




sech 
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It is useful to ask what experimental set-up could test the theory of wall tunneling in this films. Stampcil discussed 
tunneling of a wall from a defect induced on the surface of the film (the ends of the wall being held fixed by artificially 
constructed "gate" defects - this avoids the problem of edges effects in the film) . 

3. Block and Neel Walls in Wires 

The experiments of Giordano et. alM~l suggest that theoretical investigations of domain wall tunneling in wires, 
with diameter d < A, might be very useful. It should immediately be emphasised that the question of the actual 
configuration of the magnetisation in such a wire is by no means simple, and suggestions that the. magnetisation profile 
may be simply modelled as either a Bloch or Neel wall have been severely criticised by AharoniLj. This point will arise 
again later in the present paper, when we come to consider the experiments of Giordana.ei;. al. These experiments were 
analysed by the authors using the "displacement tunneling" model of wall tunneling3'El just discussed; the influence 
of both phonons and magnons on the tunneling was ignored, in conformity with the results derived by Stamp. 

Some rather more exotic tunneling processes have also been suggested for the cases where a simple Bloch or Neel 
wall configuration does occur in wires; thesej-ssfeije given the names "chiralityj-tunneling" and "Bloch coherence" . These 
possibilities were raised by Braun and LossEjO and by Takagi and TataraE2l Since we will argue later in this paper 
that chirality tunneling and Bloch coherence are highly unlikely no matter what the wall looks like, the exact nature 
of the magnetisation profile is probably not important in what follows (on the other hand, it clearly is important in 
the analysis of Giordano's experiments). 

We consider first the possibility of chirality tunneling. For a Bloch wall, chirality tunneling corresponds to a 
tunneling of the sense of rotation of the wall (from clockwise to anti-clockwise, or vice- versa). The dynamical variable 
is '/>(i), with tunneling being from (j) — ±7r/2 to (j) — with 9 constant, corresponding to the domain wall 

magnetisation rotating out of the easy plane. The tunneling splitting is very small - thus, for a wall in Ni wire, the 
tunneling splitting was estimated to be 0.1 MHz (ie., 5/iK, with a crossover to tunneling below Tg ~ 0.3 mK. 

For a Neel wall the situation is similar, except that the energy barrier is smaller, and so A,^ has been estimated to 
be as high as 80 mK. In both these cases the walls are assumed to be extremely small, containing perhaps 10^ spins 
at most. 

Braun and Loss also considered band or "Bloch" coherent motion of a 1-dimensional Neel wall (which they actually 
refer to as a Bloch wall). In the case of a Neel wall tunneling, the Berry phase of the wall enters in an important way; 
walls with integer phase have a the usual Brillouin zone, but walls with 1/2-integer phase have a halved Brillouin zone 
and the Bloch band is split into two bands, with opposite chirality. If one ignores all possible sources of decoherence 
from the environment, then a very small (10** spins) wall in YIG wire was estimated to have a Bloch bandwidth of 
Ad ^ 80 mK. It was claimed that this motion should be observable when T < Tc ~ 80 mK. 



III. COUPLING TO PHONONS AND NUCLEAR SPINS 

As noted in section I, all quantum processes are strongly influenced by the coupling to the environment. Thus the 
formulae given in section II are a rough guide only; to properly determine the tunneling rate (or whether tunneling 
occurs at all) one must first include dissipative effects. This is of course even more true of any inelastic quantum 
diffusive dynamics, which only exists because of the coupling to the environment. 

We now introduce the 2 different environments that will be studied. We begin with a few general remarks and 
useful formulae for oscillator baths and spin baths. We then go on to derive the various couplings between a domain 
wall and phonons. This is done for both 1-phonon and 2-phonon couplings. Finally, we derive the coupling between 
a domain wall and the nuclear spin bath , and show how this yields a slowly fluctuating spatially random potential 
acting on the wall. 

A. Oscillator Bath and Spin Bath Environments 

As already remarked in section I, we will be ignoring the effects of magnons, photons and electrons in what follows; 
in all cases this is because their dissipative effects are weak. We are nevertheless left with 2 important environmental 
couplings. Before giving a detailed treatment of these, we make some general remarks and recall some general formulae 
that will be relevant in the later technical discussions. _ 

The phonpn bath is representative of the oscillator bath model discussed by Feynman and VernonO, and Caldeira 
and LeggettES. A subtlety that will arise here is that the individual oscillators do not necessarily represent individual 
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phonons: they can also represent pairs of phonons (an analogous situation was encoiiHtered in the discussion of 
magnon coupling to domain walls, where the oscillators represented triplets of magnonsdu). 

The nuclear spins on the other hand, are a good example of a "spin bath" , in which each environmental mode behaves 
as a 2-level systenY(wt|h associated 2-level Hilbert space) . In general the spin bath can behave quite differently from 
any oscillator batbu3E3. 

The interaction between the nuclear spins and the phonons is utterly negligible. Thus the 2 baths act separately 
on the domain wall. Their effects are quite different and depend on the different nature of the 2 baths as follows: 



(i) Oscillator Baths : This model assumes the environment can be described as a set of non-interacting harmonic 
oscillators, with each environmental mode weakly coupled to the system. For N delocalised environmental modes, the 
coupling between each oscillator and the system ~ 0{1/^/N) and the interaction between oscillators is ~ 0{1/N). 
The coupling between system and environment is linear in the bath coordinates. We emphasize that the connection 
between the oscillator modes and the modes that might appear in some microscopic model of the system can be 
non-trivial and non-linear. ProvidecLsuch an oscillator bath exists, one may write a general Lagrangian now known 
as the Caldeira-Leggett Lagrangiantj, 

LcL = yQ' + V^(Q) + + +^kXk') + 5I[^k(Q, Q)xk + Gk(Q, Q)xk] + <i>(Q, Q) (3.1) 

k k 

where the environment is represented by the set {xk} and <i>(Q, Q) is a counterterm used when a phenomenological 
equation of motion is available; it simply cancels the shift in the potential introduced when the system is coupled to 
the oscillators. 

An important simplification occur if an expansion of i^k(Q, Q) and Gk(Q, Q) to first order in Q and Q is possible. A 
further set of transformations on the batho eliminates the coupling of the bath to Q, and then writing Jk(Q) ~ GkQ, 
the Lagrangian is written as 

Lo + iB + jQ^CkXk-iq^^-^ (3.2) 



Lc-L — 

z - — ■ z 

k k 



where Lq is the system Lagrangjaa and Lb the oscillator bath Lagrangian. The reduced density matrix of the system 
then defines the effective actiont^O Sef f = Sq + ASef / , where Sq is the system action andllj 

ASeff dr I dr'a{r - r')(Q(r) - Q(r'))^ (3.3) 



with the kernel 

rv(T — t'] = 

27r 

where the environmental mode propagator D[uj,t) and spectral function J{lo) are 



a{T-T') = ^ I dujJ{Lo)D{uj,\T-T'\) (3.4) 



D{uj, t) = cosech(tj/2r) cosh(w(^ - |r - t'|)) (3.5) 



^M-?E^'^(^-^k) (3.6) 
2 ^ mk^k 

k 



Alternatively, one can write 



1 iuj„{T-T') 
-r') = -rTyj dujLuJioj) ^ , , (3.7) 



where cun = ^irnT is the bosonic Matsubara frequencyEl. 

It should be emphasised that these effective Lagrangians, and the functions a(T) and J{i^), refer to a Hilbert space 
for both system and environment that has already been truncated to low energies. High energy environmental modes, 
as well as higher energy states of the system itself, are incorporated into the Lagrangian in the form of renormalised 
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couplings. Consequently the high-energy form of J{ijj) typically has some smooth cut-off (for example of the form 
J(w) ^ exp(— w/wc))- At low frequency, J(w) often has power law form. Most dissipation is caused by any term in 
J{ijj) of Ohmic linear form, ie., for J{uj) = rjoj where is the classical friction coefficient. In the absence of Ohmic 
dissipation, there will always be supcrOhmic terms of the form J{uj) ^ uj'', with s > 1. Occasionally, one may 
also have to worry about a sharp jump in J{uj) at some finite frequency, usually imposed by a gap in some set of 
environmental modes. However in this case (a) other modes will always give contributions to J{lu) for uj below the 
gap energy, and (b) higher-order couplings to. the same gapped modes are typically ungapped. Thus, for example, in 
the case of magnons coupled to a Bloch wall3, even though the magnon spectrum is gapped, one finds that coupling 
to pair of magnons gives an ungapped superOhmic ( J(w) ^ i^^) term, and coupling to triplets of magnons gives an 
ungapped Ohmic term. Note, however, that such non-linear couplings to multiplets of bosonic excitation are always 
temperature-dependent, ie., J(a;) J{uj,T). 

In a tunneling situation, the strength of thcj-eoupling to the environment is determined by the dimcnsionless 
parameter at — ri/2Myj^Q for Ohmic couplingij, and we can define /?t = P^^/Myj in the case of superOhmic 
dissipation with s = 3, when the spectral function J{uj) — Puj^. 

In the presence of friction, the crossover temperature from quantum tunneling to thermally activated-relaxation is 
decreased. The new crossover temperature can be estimated in the case where J{uj) is T-independentJI3o 

^ Tolil + a^t^/^ - at] (3.8) 
for Ohmic dissipation. For superOhmic dissipation, in the weak dissipation regime (/3t <^ 1), Tc — To(l — /3t/2), while 

1 /3 

for strong coupling Tc — Tq/A (again for a T-independent J{uj)). 

(ii) Spin Baths: In a wide variety of cases the environmental modes each have a finite Hilbert space. In the simplest 
case where this Hilbert space is 2-dimensional, each environmental mode is equivalent to some 2-level system, or to 
a spin- 1/2. Examples of this include nuclear spins, paramagnetic impurities, as well as a variety of more subtle 
2-level systems existing in glasses or disordered solids, often associated with defects. One should also note the various 
attempts tp_model the behaviour of dissipative quantum systems in terms of a set of 2-level "Landau-Zener" degrees 
of freedonxj. 

Just as for the oscillator bath environment in Eq. (^), one may write a a general description of a system interacting 
with a spin bath. It is more convenient to use a Hamiltonian formulation of this; the most general Hamiltonian form 
is then 

HsB ^ Ho{P,Q)+HBi{(rk}) + H,nt{P,Q,(Tk}) (3.9) 
where P is conjugate to the system coordinate Q, Hq describes the system, and the spin bath has an Hamiltonian 

HBiM) = E • '^fc + ^ E E Kk'^<'>^ (3.10) 

in which the (T^ are Pauli spin matrix dynamical variables (with k = 1,2,... N), and the h^'s are "external fields" 
which may exist. The couplings V^jS between the spins are often assumed to be much weaker than the characteristic 
energy scales of the system (for nuclear spins they are the very weak internuclear dipolar interactions; \Vkk' \ ~ 10~^ K 
or less). 

The interaction Hamiltonian has the form 

i?„,,(p,Q,^fc}) = E (^fe(P'Q)^fc + ^(^fc^(P'QK + n"(P'QK)) (3-11) 

k ^ ^ 

containing both "diagonal" couplings -F'^(P, Q) which polarise the environmental spins) and "non-diagonal" couplings 
F^±(P, Q) (which flips them). 

Whereas in the case of oscillator baths, averaging over the oscillators is conveniently done using path integrals 
(to produce an "influence functional" ) , the lack of a well-defined classical path for the 2-level system makes such an 
approach infeasible here. Moreover, we stress that the couplings -^^'(Q) are in general not small. This is in sharp 
contrast to the oscillator bath model, where the couplings Fk(Q) ~ 0{N~^l'^')\ in the spin bath the couplings are 
independent of N (so that their effects become more and more serious as N increases). For hyperfine interactions, 
F^(Q) can be as high as 10 GHz (~ 0.5 K), for Ho nuclei - thus one can not always assume that a perturbative 
treatment of the -F'^(Q) is valid. 
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To understand the effect of the spin bath on the dynamics of Q, we make the assumption that F^{P, Q) F^(Q). 
This is true for all cases studied so far (including domain walls). In this case, we may distinguish 2 effects. 

On the one hand, the diagonal couphng U{Q) = F|(Q)cr| acts as a random potential in the coordinate Q; we 
shall see this in detail for the case of a single domain wall. It is useful to define a density of states W{U) for the 
potential U. In an ensemble of potentials (or of different systems), and assuming a typical value ojq for F|(Q), W{U) 
is Gaussian in shape, of form 

W{U) - {2ttE^)-^/'^ exp[-[/V2£^o] (3-12) 

where Eq ^ N^^'^loq is the gaussian half-width. In general the potential will fluctuate in time, because of transverse 
spin relaxation, at a rate governed by the nuclear Tj""'^ (we will ignore Ti processes in this paper, assuming that 
Ti is very long). Thus i7(Q) also varies in time, and fluctuates, for any particular Q, throughout the domain of 
W{U). If J7(Q) were static, it would simply block the quantum dynamics of the system, by removing the degeneracy 
between initial and final tunneling states (we call this effect "degeneracy blocking"). Its fluctuations in time mean 
that "resonance windows" are occasionally opened between such states. Notice that the effect of ?7(Q) is entirely 
elastic - no dissipation is involved. 

On the other hand the non-diagonal couplings F^(P, Q) cause environmental spins to flip under the action of the 
system. Thus, whereas the main effect of F^(Q) is to completely alter the effective potential acting on the system, 
the main effect of the non-diagonal couplings is to cause transitions in the spin bath, leading to dissipation and/or 
decoherence. Unlike the oscillator bath, dissipation and decoherence are not necessarily linked in the spin bath - 
indeed it is possible in certain cases to have decoherence without dissipation. Because the main effect of transitions in 
the spin bath is the absorption of a random topological Berry phase by the bath (thereby randomising the quantum 
phase associated with motion in Q), this decoherence is called "topological decoherence". 

In this paper we shall develop the theory for a domain wall coupled to nuclear spins, but we shall ignore the 
dissipative effects coming from nuclear transitions. While these are not necessarily small, the "degeneracy blocking" 
effects, coming from the diagonal coupling ^'^(Q) are often considerably larger, particularly for weak hyperfine 
potentials (such as in Ni or Fe materials). Their main effect is the generation of a random potential acting on the 
wall, whose effects must be understood before dissipation is considered. 



B. Phonon Couplings to the Wall 

As noted above, it is not sufficient to consider only 1-phonon couplings to the wall; we must also consider non-linear 
2-phonon couplings. These will be analysed using standard continuum elasticity [tiieory for acoustic phonons coupled 
to the magnetisation; this rather intricate subject is reviewed by de LacheisserieEa 

The Euclidean Lagrangian of the acoustic phonons in a material of mass density is 

£p{x,T) = ^PvUi + ^CijkiUijUki (3.13) 

where Ui — Vi — represents the phonon's field, ie., the displacement of the atoms with respect to their equilibrium 
position, and Uki = {dkUi + diUk)/'2 is the strain tensor. The potential energy of the field is given by the elastic tensor 
Cijki, with typical values C ~ lO^^J/m . The indices refer to the directions of the displacements, with the summation 
convention over repeated indices applied. We begin by considering an isotropic elastic energy, which is the simplest 
form 

Cijki = KSijSki + fJ-e{5ikSji + SiiSjk) (3-14) 

where fj,e and Ae are the Lame constants. The tensors that are going to be considered in this paper are all symmetric 
under the exchange of two indices forming a pair {ij} and it is convenient to use the abbreviated notation 

{11} ==1 {22} = 2 {33} = 3 

(3.15) 

{23} = 4 {13} = 5 {12} = 6 

We will thus represent any pair of indices by a letter a — {aia2} = {a2ai} = 1, ...,6. This simplifies the notation 
tremendously. 

We can also consider the phonons of a cubic structure. In this case, the non-zero elastic constants are Cn = C22 = 
C33, C12 = Ci3 = C23 = . . . and C44 = C55 = Cee. It reduces to the isotropic case if Cn — C12 — 2C44 — 0, in which 
case C12 — Ae and C44 — 2/ie. 



11 



It is observed experimentally that the shape of a fully magnetised ferromagnetic substance is different from the 
shape in the unmagnetised state. This phenomenon, called magnetostriction, can be related to the strain dependence 
of the anisotropy energy. There is thus an interaction between the direction of the magnetisation and the strain 
tensor of a solid. In many cases, this interaction is linear in the strain tensor. The simplest type of interaction that 
respects time-reversal symmetry is thus of the form Uijifikrhi, with j, j, k and / arbitrary directions. The interaction 
between phoapns and the magnetisation is mediated by the first order magnetoelastic tensor, Aijki , as discussed in de 
LacheisserieES. It is also possible to consider an interaction that is quadratic in the strain and inithe magnetisation, 
which is described [by the second order magnetoelastic tensor Rijkimn first introduced by Masoro and also discussed 
by de Lacheisseriala. The general form of the interaction between the phonons and the domain wall is given by 
integrating over these local tensor couplings, throughout the sample, for a given wall profile, and then subtracting off 
the result that would have been obtained in the same integration in the absence of the wall. Thus, if the "vacuum" 
magnetisation profile is written as m*^, the general form of the interaction Lagrangian will be 



Lint = - d rAijkiUij{r){mk{r)mi{r) - m^,mi) + RijkimnUij{r)Ukiir){m„,{r)mn{r) 



(3.16) 



up to 2nd order in the phonon strain variables. The "vacuum" (no wall) term is just a constant (the integrated bulk 
magnetoelastic stress energy for the entire sample), which we subtr act o ff in all our calculations. 

^From now on, we will use the definitions of the indices in Eq. ( 3.15| ) to discuss this Lagrangian. The first order 



tensor has been extensively studied in the context of magnetostriction. The second order tensor is encountered less 
frequently, but has been used in connection with the so-called "morphic effect", ie., the change in sound velocity as 
a function of the direction of the applied fieldEj. Both of these couplings come from an expansion of the exchange 
energy in terms of the displacement of the atoms. In our context, A and R then correspond to 1- and 2-phonon 
interaction terms respectively. Note that 2-phonon ms^us that the interaction involves two phonons simultaneously, 
a process quite different from two 1-phonon processescZlE2l. 

For a crystal with orthorhombic symmetry, there are 12 non-zero matrix elements: An, A22, A33, A12, A21, A13, 
A31, A23, A32, Aqq, A55 and A44. If the symmetry of the crystal is reduced to cubic, then An — A22 = ^33, 
A44 = A55 = Aqq and A12 — A21 — A31.... The second order magnetoelastic tensor is obviously very complicated, 
but we will only use i?iii — R222 — R333 with all other Rabc equal to zero. This approximation is justified by 
the fact that these three components are generally of the same order of magnitude, while being at least two orders 
of magnitude IsH^fe^ than the other cop^cients. Several measurements and calculations of these coefficients hq 



been performecEl'Eil. In YIG— these areEJ: Aab ^ 10 J/m 
CT ~ 3 X lO^m/sec. In nickelEa: Aab 10* J/m^, Rm 
We use the Fourier transform of the phonon field 



i?iii ~ 10^ J/m^, with the transverse sound velocity 
10^° J/m'^, and ct ^ lO'^m/sec. 



u(r,T) 



n q 



(3.17) 



to write the action corresponding to the interaction Lagrangian, Eq. (3.16) as S'/[u, rh] = 5*1' ■* [u, rh] + Sf^' [u, m] with 
5f [u,m] = ^AabTY^Y. / dTe'^-'i^^'>e-'^"^Mb{-ci) Ki^a, (q, *^„) + 9a.Uax (q, ic^„)] (3.18) 
representing the action coming from 1-phonon processes and 

r-l/T 



:(«'0 r 



Sf'\u,rh] = -^RabcT^Y^Y. J dTe*('^+*'')-QMe-^(""+'^"')^A^,(-k-k') 

nn' kk' ^ 

[ka-,Ua2{K i^n) + ka^Ua^Ck, iuJn)] [k'f^^Ub^ik' , iuj'„) + fc^^ Mfc^ (k' , iw'j] 



(3.19) 



corresponding to 2-phonon processes. Q(t) is the position of the wall and the summation convention over repeated 
magnetoelastic indices a = {0102} is in effect. 

In these two expressions, the profile of the magnetisation is included in the magnetisation form factor A4a defined 

as 



(3.20) 



The magnetoelasticity of thin films differs considerably from what is observed in the bulk@. The elastic properties 
and the symmetry of a film will in general be different than those in bulk, and the strain associated with the surface 
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will also play a major role. The magnetostriction constants in general are dependent on the thickness of the film and, 
furthermore, in many cases when discussing the simple magnetostriction, it is impossible to consider only the linear 
magnetostrifitkttt. It turns out that it is strain dependent and that the second order magnetoelastic constants play 
a major roldlSCJ. It is obvious that a proper discussion of the magnetoelastic dissipation requires a v ery d etailed 
knowledge of the material being used. To keep the discussion general, we still use the Lagrangian Eq. (3.16) where 
now the phonons will be 2-dimensional, with a surface density ps and the elastic and magnetoelastic constants have 
units of J/m^. Similarly, in 1 dimension, we consider a linear density pi and the constants have units of J/m. More 
detailed models could be devised if necessary. 

Finally, a special comment must be made about the terms in the magnetisation profile appearing due to the finite 
velocity of the wall, ie., coming from the demagnetisation energy. They represent a coupling to the environment 
that is proportional to the velocity of the macroscopic coordinate. The standard way to deal with this coupling is 
to introduce a total tinae, derivative term in the action and to perform a canonical transformation to a new set of 
oscillators (c.f. LeggettE3). In our case, it does not appear that this can be done easily since the coupling is not 
"strictly linear", but is a complicated function of Q. We will first build the effective action of the magnetisation 
including the velocity term. This will allow us to see in which cases exactly this coupling to the velocity is important. 
Once its physical relevance or irrelevance is established, it is easy to go back to the original Lagrangian and to deal 
with this term directly. 



C. Nuclear Spin Coupling to the Wall 



Nuclear hyperfine effects vary enormously between magnetic systems. The weakest is in Ni, where only 1% of the 
nuclei have spins, and the hyperfine coupling is only ujq — 28.35 MHz (~ 1.4 mK). On the other hand, in the case of 
rare earths, wq varies from 1 to 10 GHz (0.05 to 0.5 K). In this latter case the hyperfine coupling energy to a single 
nucleus may be comparable to the other energy scales in the problem ! 

To derive an effective interaction Hamiltonian, we consider our system of ferromagnetically ordered spins to coupled 
locally to N nuclear spins \k at positions {k = 1, 2, 3, ... N), which for a set of dilute nuclear spins (where only one 
isotope has a nuclear spin) will be random. The total Hamiltonian for the coupled system is then 



N 

E 

fc=i 



t^fcSfc ■ Ife 



2 E E ^kk'^k^k' 



(3.21) 



where Hm is the electronic Hamiltonian for the magnetisation (Eq. (2.1)), written in terms of the electronic spins 
at the sites where there happen to be nuclear spins, and oJk is the hyperfine coupling at r^; V^^^, is the internuclear 
dipolar interaction,, with strength \Vl^^\ ^ 1 — 100 kHz (0.05 — Q.hpK). In terms of the continuum magnetisation 
M(r), we have 



H = Hr„ + Y,LOk I — <5(r-r,0[M,(r)/| + (Af.(r)4-+M,(r)4^)] + i^^F,f4"/f, (3.22) 

k=l ^3 ^ k k' 

The resemblance between thi s form and the spin bath coupling form given above is evident - however, it is often 
more convenient to rewrite Eq. ( 3.22 ) by separating out the slowly- varying wall profile Mvi/(r) from the "fast" magnon 
fluctuations which ride on top of this. Thus we write 



M(r) =Mw(r,Q) + JM(r) 



(3.23) 



where MvK(r, Q) is the wall profile for a wall centered at Q (given in Eq. (O) for a Bloch waU and Eq. ( p.lTD for 
a Neel wall). We then introduce local axes {u,v,w), with ic^, parallel to_Mvi/(r, Q), and Xu, perpendicular to it, 
and to each other (these axes were called xi,X2,X3 in previous paperaafl). Then we quantise the ^M(r) in terms of 
local magnon operators 5(r) and b~^{r), as usual: 



5A4(r) = -47g6+(r)fe(r) 



SM+ir) = (47,Mo)i/2 (^1 - ^6+(r)&(r)) 6(r) 
<5Af_(r) = (47,Afo)^/'&+(r) (l - ^fe+(r)6(r)' 



(3.24) 
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where 5M±(r) = (5M„(r) ± SMy{r). The interaction term in Eq. (3.22) is now written for a single domain wall by 
again subtracting off the "vacuum" magnetisation profile Mq obtaining in the absence of the wall, to get a sum of 
longitudinal and transverse couplings: 

H^ntiQ, {Ik}) =ll^r^nt +^ H,r.t (3.25) 

N „ ,3 

lli/^nt = / —Kr ~ rfe) ((A%(r, Q)I^ - Mo • Ife) + I]:SMUr)) (3.26) 



k=l 



1 ^ f (P 

^H„u ^T^Y.^'k —S{r - Vk) {5M+{v)I- + <5M_(r)/+) (3.27) 

^ fc=l ^3 

To get an idea of the size and effect of each of these terms, let us choose as example a Bloch wall, supposing for 
simplicity that ujk = wq for each nuclear spins (in reality the {ujk\ will be spread out by the internuclear V^'^f , as well 
as by transfer hyperfine couplings to non-magnetic ions). 



Supposing also that the Bloch wall is static, we have, from Eq. (2.4), for a wall centered at — Q, a, static 
longitudinal interaction 

X3- Q 



(Q, {I.}) - ^ E / d'rSir - r,) ((^1 - Ctanh 



A 



B 



5m.u,{r) (3.28) 



and a transverse interaction 

N 

-L rjBloch 



Hg[-\Q, {I.}) = ^ E / d'rSir r,) (xscch (/+ + /") 

+ {6m^{r)I+ + 6m+{r)I^) ) (3.29) 



where Sm = SM/Mq. 

We can extract from Eq. ( 3.28 ) and ( 3.29| ) the following 3 terms: 



(i) Terms of the form I^Smp couple nuclear spins to magnons; this interaction is considerably enhanced inside the 
wall where the magnons are gapless, and is known to be important in the discussion of the NMR relaxation rate inside 
the wall. 

(ii) We also have a longitudinal tanh[(a;3 — Q)/Ab] term in "-ffmt which is diagonal in the local nuclear spin basis 
states (with axis of quantisation ic^ along Mvy(r, Q)). This term gives us a random potential when summed over all 
spins, which we write as 

UiQ) = ^ E/ d'r6{r - r,) (1 - Ctanh (^)) (3-30) 

Suppose we use the "high-T" limit, where kT ujq, so that the expectation value (I^) is zero in any direction. 
Since the nuclear spins are completely uncorrelated (VJ^^f is negligible), a volume containing TV nuclear spins will have 
root mean square polarisation ~ N^^^. It then immediately follows that for a system in which a fraction x of all the 
states are occupied by nuclear spins Ifc at sites fc, the ensemble averaged correlation between U (Q) for 2 different wall 
position is 

CuuiQi - Q2) = ((f/(Qi) - U{Q2)f) - u;Is^I^ANq,_q, 

^^^M^] alxS^{Q,~Q2) (T»c^o) (3.31) 

where s = |sfe|, / = |Ifc|, S'^, is the area of the wall, and AN{Qi — Q2) ^ {xSi^/aDlQi — Q2) is the mean number 
of nuclear spins in the volume swept out by the wall in going from Qi to Q2- Thus U{Q) exhibits "random walk" 
behaviour over the sampleH, as represented on Fig. (|[). 
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Eq. (3.31) is for the high-T limit. In the case of very low T, when T <C Wq, the nuclear spins will tend to line up in 



the field due to the electronic moments. In the low-T limit, they will all be aligned according to the wall configuration, 



and Eq. (3.31) is replaced by 



Cuu{Qi - Q2) - alx'SliQi - Q^f (3.32) 

Thus the wall is trapped in a potential which increases linearly (on average) in both directions away from the wall 
centre. In reality, both results ( |3.31 ) and ( |3.32 ) are only valid for times <^ T2 ; for larger times the internuclear dipolar 



transitions cause random fluctuations of the hyperfine potential. 

(iii) Finally, we have the sech[(a;3 — Q)/Xb\ term in Eq. ( 3.29| ), acting only inside the wall; this simply tells us that 



in the chirally- rotated frame oiTs/lwi'^, Q), the nuclear spins quantised along the "vacuum" Mq are no longer diagonal. 
Because of this term, a domain moving through a field of nuclear spins can fiip some of them as it passes, thereto 
causing both topological decoherence and dissipation in general. From the solution of the "central spin" problerrO 
one may estimate that the mean number A/ of flipped spins caused by a single excursion between 2 points Qi and Q2 
(sweeping out AA^ nuclear spins) will be given by 

where is the typical frequency involved in the motion from Qi to Q2 {^a being the bounce frequency for tunneling 
motion), provided wq ^ f^o- 

To get a feeling for the numbers, consider the "high-T" limit fc^T ^ wq- Now imagine a wall of area (300 x 300)A^, 
which sweeps out a length Qi — Q2 ^ 300A, and a volume containing 10^ electronic spins (ie., qq ~ 3A); suppose 
also that 5 = |Sfe| = 1 and I — \lk\ — 1/2. We may imagine 2 extreme cases: 

(a) Suppose uiq is very small, eg., let ujq = 1.4 mK, and let x = 0.01 (only 1% of nuclei have spins); these are values 
appropriate to Ni. Then the Gaussian mean change in U{Q) for the wall will be \AU\ ^ 0.3 K, in the high-T limit 
(ie., T ^ lAmK; all experiments up to now have been in this limit). If we assume that Hq ^ 2 x 10^" Hz (typical of 
wall tunneling), then Xj ~ 0.025, ie., in most cases no nuclei will be flipped at all. 

(b) Suppose now loq — 0.2 K, a value more typical of some rare earth magnets, and that x = 1. We than have, 
in the high-T limit (now for T > 0.2i^), that \AU\ 400 K, and A/ ~ 10^ (each nucleus has a probabihty ~ 0.1 of 
being flipped) ! Even if we reduce the displacement in this case from 300A to only 3A (ie., the wall makes only a 
single lattice displacement), we still find \ AU\ ~ 40K and A/ ~ lO'^. 

Notice furthermore that the low-T limit fc^T <ti luq, where all nuclear spins line up with the local magnetisation, 
is easily attainable for rare earth systems, in current experiments. In this case the effects are spectacular; for the 
present example, with luq = 0.2 K, a displacement of 300A leads to an energy change |AJ7| ~ 20 eV!! Thus the wall is 
completely locked to the nuclear polarisation when fc^T <C 0.2 K. 

Thus in the first case, of extremely weak coupling, we see that dynamic effects in the bath can be ignored - the sole 
function of the nuclear bath is to provide a weak random static potential. On the other hand in the case of rare-earth 
magnets, with very strong hyperfine coupling, the nuclear effects are quite enormous, and moreover it is crucial to 
take account of nuclear transitions in analysing the wall motion. 

In this paper we will concentrate on the weak coupling case, because it is simpler and more relevant to the Ni wire 
experiments of Giordano et. al. We shall deal in section V with the residual effects of the internuclear coupling, which 
drives the slow (on the scale of rig ^) fluctuations in time of U{Q). 



IV. EFFECTIVE ACTION OF A DOMAIN WALL COUPLED TO PHONONS 



In this section_we obtain the effective action of a domain wall interacting with the phonons through the magnetoe- 
lastic interactionL3; we completely ignore the nuclear spins in this section. We consider both 1-phonon and 2-phonon 
processes, asjdescribed in the last section, and we relate their effect to the phenomenological oscillator bath of Caldcira 
and LeggettEj. As expected 1-phonon processes give rise to superOhmic dissipation with J(w) ^ uj^. We find however 
that the dominant contribution to the dissipation comes from the coupling between the velocity of the wall and the 
phonons. 

The second order magnetoelastic tensor brings in 2-phonon mediated Ohmic dissipation, with a temperature de- 
pendent friction coefficient t]{T) ^ (-^vheije-id is the effective dimensionality of the phonons), just as with the 
quantum diffusion of a particle in an insulatoiEd. In this case however, the coupling to the velocity of the wall is 
unimportant. As part of our the discussion of 2-dimensional 2-phonon processes, we discuss the general case of the 
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dissipation of a Neel wall. This represents a straightforward extension of the analysis of a Bloch wall. The same is 
true for the dissipation in a wire, where the phonons are effectively 1-dimensional. 



A. Effective Action of a Bloch Wall 



We are interested in integrating out the phonons so as to obtain an effective action for the magnetisation. The 
effective action is defined as 



-So[m 



-So[u]-S7 [m,u] 



e-So[rh],g 



jD[u]e-SoH 

-Sj[m,u]\ 



(4.1) 



where 5*0 is the action of the phonons and Si is the magnetoelastic interaction , as described in the Introduction. We 



recall that the inter action action Si = s'f'^ + Sy , representing 1- and 2-phonon processes, is expressed as the sum of 



Eq. (|3.18D and Eq. ( |3.19| ). The uniform magnetisation, in absence of a domain wall, is in the xi direction. The form 
factor of the magnetisation, appearing in the action Si is thus 



(Tre 



''(TOai(x)TOa2(r) - Sa^.lSa^,!) = (27r)^(S(gi )(5(g2 ) A^a (^s) 



(4.2) 



This last form coming from the one-dimensional properties of a Bloch wall. The magnetisation depends only on 
the coordinates normal to the wall and this restricts the momentum in the plane of t he wall to be zero. Using the 
description of a Bloch wall in term of the component of the magnetisation Eq. (2.4), it is easy to compute these 
different terms. Up to second order in the velocity of the wall, they are: 



■^f (93) = -7r(j3A^cosech(7rq3A/2) 



(4.3) 



M^{q:i,T) = ( 1 - ) 7rg3A2cosech(^(j3A/2) (4.4) 



^0 



M^{q3,T) = :2-^7rg3A^cosech(^(73A/2) (4.5) 
4Cq 



Mf{q3,T) = ^^q3A2cosech(7rg3A/2) (4.6) 
2co 



Miiqs, r) = -^^nq3Xhec^nq3X/2) (4.7) 

ZCq 



Miiqs, t) = -iIi- j 7rq3A2sech(^q3A/2) (4.8) 



where A = {J/K^^y^^ is the domain wall width and cq is the Walker velocity, defined in Eq. (2.5). There are 
essentially two types of coupling; those proportional to cosech(7r(73A/2), which go to the constant 2A as 93 — + 0, while 
those proportional to sech(7rg3A/2) ^ as 53 ^ 0. This already indicates that dissipative terms involving Aii, M2, 
Ms and AI4 will be more important than those coming from A4q and Ai^. 

It should be noticed that there is a natural cutoff momentum kc = 2/(7rA) coming directly from the wall structure. 
This can be converted to a cutoff frequency ojc = 2ct/{ttX). For a wall of width A = 500A, and with a transverse 
sound velocity ct ~ lO^m/s, this gives a cutoff frequency lUc = 10^^ Hz. Walls with a very small width wil l hav e an 
even larger uj^. This cutoff frequency is usually much larger than the bounce frequency r^o; defined in Eq. ( 2.15| ), at 



which the bath must be cutoff from the truncation procedure. This makes it possible to use cosech(w/ajc) ^ '-^c/'-^ 
and sech(u;/ijjc) ^ 1- 
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The exact form of the magnetisation profile is actually not too critical; we are mostly interested in the behaviour 
of these form factors at \ow-q. They are essentially an integral of the components of the non-uniform magnetisation 
out of the easy axis. The spatial variation of these terms is actually what defines the width of the wall, so it is quite 
clear that the relevant form factors for the dissipative dynamics of the wall behave as 

\im M{q)^BxX (4.9) 

where B\ is a coefficient of 0(1) depending on the precise profile of the magnetisation and A is the length scale over 
which the the magnetisation changes from one stable configuration to another. This should be accurate provided that 
the curvature of the wall is weak, of course. _ 
The two-point function of the phonon field is decomposed as usual into a transverse and longitudinal partH 

{Ui{p,iuJn)Uj{-p, -iuJn)) = {Sij - PiPj)G'^ {p, iuJn) + PiPjG^ {p, iUn) (4.10) 

with 

G^'^(p,zc.„) = -^ 1—^ (4.11) 

where — He/ Pv and c| = (Ag + 2iie)/ Pv arc the transverse and longitudinal sound velocities. 

The effective action is found from a simple cumulant expansion. The first term, {Si) is ignored: the 1-phonon 
term {S^'^) is zero due to (u) = 0, and the 2-phonon term is independent of Q{t), simply adding a contribution to 
the potential energy of the phonons (morphic effect). The ratio of the second-order magnetoelastic constants to the 
elastic energy is generally less than 1% and we can simply ignore this renormalisation in the sound velocity. The first 
dissipative contribution to the effective action is thus {Sj)/2 = {{SY^^Y) /'^ + ((*S'j^^)^)/2 which gives 

A5e// = l-AabAedTy^Y] /'^^dTdT'e'-"(^~^')e''?-«5(^)-«(^'»Alf,(-q)Ald(q)gac(q,it^n) 

n q ^ 

+ \Rai>cRdefT^ J2 E dTdT'e^(""+""')(^-^')e^('=^+'=3)-(Q(-)-Q(-'))x^(_k - k') 

kk' "'0 



nn' 



xMfik + k') [Gad (k, iuJn ) Gbe (k' , iuJn') + GaeiK iUn ) Gbd (k' , iw„0] (4-12) 



where Gab = |(fcai ^61^0262 (k, *w„) -I- Perm). The index a = (ai, 02) is the pair of the indices introduced in Eq. ( 3.15| ) 
and Perm indicates the permutations ai <-> 02 and 61 <-!■ 62. A third term, of the form GabGdei is generated in the 
effective action, but it does not couple to the position of the wall, simply corresponding to a further renormalisation 
of the sound velocity. Since it is not a dissipative term we do not include it in the effective action. It is now possible 
to examine the different contributions to this action. 

B. 1-Phonon Contributions to Action 

The 1-phonon part of the action is very simple. Summing over the frequency, and separating the action into a 
transverse and a longitudinal part, we obtain 

ASlP'} = ASlf/ + ASlP'}'' (4.13) 



CT Pv Jq 27r Jo 

[Al^M5{-q3)M5iq3) + Al^Mi{-q3)Mi{q3)] (4.14) 



ASlff = ^^l ^1 dTdr'cos{q^{Qir)~Qir')))q,DicTq3Ar~r'\) 



A qlph,L 



'icT Pv Jo 27r Jo 

A3aA3bMai-q3)Mb{q3) 



drdr' cos(<Z3(Q(t) - Q(t'))) 93 D{cLq3, \r ~ t' 



(4.15) 



17 



where D{uj,t) is the boson propagator, defined in Eq. (3.5). This simple form of the effective action is obtained by 
taking advantage of the fact that for a Bloch wall gi = 92 = 0. Next, inserting the expressions for the Aia's, we 
obtain terms that contain only the position of the wall, of the form cosw((5(r) — Q{t')). However, there will also 
appear terms like Q^(t) cosuj{Q{t) — Q{t')) and Q{t)Q{t') coslj{Q{t) — Q{t')) coming from the velocity-dependent 
terms in Eq. (U) to (U). We can rearrange Q(T)g(r') as (Q^(r) + Q^t) - {Q{t) - Q(r'))^)/2 which brings a 
constant shift in the domain mass and a dissipative term, which dep ends on the velocity of the wall. We then obtain 
an effective action that can be compared immediately with Eq. ( p.Sj ), coming from the Caldeira-Lcggett modell3 : 



ASeff[Q,Q]^l [ / dT'{QiT)-QiT')faL{r^T') 

^ Jo Jo 

.1/T .1/T jTj .1/T 

dr dT'{Q{r)~Q{T'jfaT{T-T')-— Q\t) 
a Jo ^ Jo 



(4.16) 



where 



aL = ^-^(^32 - Asif [ duJLu^cosech^ f?^) k - ^'D 
4 P c£ Jo V 2cL / 



(4.17) 



aT 



TT S X"^ 2 



dojtj cosech" 



/ TTLUX 



D{u;,\r-r'\) 



(4.18) 



M = 



3 Pv Cq 



-((A32 - ^31)^33 



(^32-^31)72) + — 



(4.19) 



This effective action is not completely equivalent to the Caldeira-Leggett action, due to the terms in Q. This is of 
course a perfectly legitimate action, but its analysis is not simple. It however allows an easy identification of the 
dominant contribu tion to the dissipation. The normal dissipative term is of the superOhmic form (^(ip4 ~ uj^) as can 
be seen from Eq. ( 4.17| ). This is what would be expected from one-dimensional one-phonon processesO. It must also 
be expected in the case of dissi pati on for a Bloch wall, since the wall restricts the phonon momentum in the plane 
of the wall to be zero (see Eq. (4.2). The dissipative term coming from the coupling of the phonons to the velocity 
is described by a spectral function linear in w, and looks at first sight to be of the Ohmic form. However, since it 
depends on the velocity of the wall, its effect is similar to superOhmic dissipation. 

In the transverse kernel arir — t'), there appears the Walker velocity cq, defined in Eq. (^.£|), and the transverse 
sound velocity whereas only the longitudinal sound velocity appears in ai(T — r'). In general, Cq < Ct < cl so 
that depending on the value of the coefficient A44, the term coming from the velocity of the wall might be the most 
important. Secondly, for a crystal of cubic symmetry, ^32 = A31 and dissipation comes only from the coupling 
between the phonons and the velocity of the wall. This can be traced back to our neglect of the influence of the 
coupling to the lattice in the determination of the shape of the wall. Taking it into account would give dissipation 
with a coefficient that would be smaller than the coefficient of the velocity-coupled dissipation by at least a factor of 
A-ab/Ccd, where Ccd are the elastic constants. For nickel, this is a reduction of ~ 10~^, while for YIG it is ~ 10"'*, 
and it can be safely forgotten in both cases. This shows again that the coupling to the velocity of the wall is quite 
important. Let us then go back and perform the analysis by considering only the coupling between the phonons and 
the velocity of the wall. 



C. 1-Phonon Coupling to the Wall Velocity 



We now look exclusively at a wall with a cubic crystalline structure or with cq < ct such that the dominant 
mechanism for dissipation is from the coupling of the phonons to the velocity of the wall. As the wall is a one 
dimensional structure (in a flat wall approximation), one has the restriction fci = fc2 = 0. The only relevant terms in 
the magnetoelastic Lagrangian are thus 

=^e*'=^^(")[A33f/3(fc3,T)A^3(-fc3)+A3iC/3(fc3,T)(Xi(-fc3)+A42(-fc3)) 

+ AMU4{h,T)MA{-k3) + A55U5{h,T)M5{-k3)] (4.20) 
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and A4i + M2 are proportional to Cf' while A^4 and A^s are proportional to Q. We are considering exclusively 
tunneling of the wall so that the excursions in position can be restricted to second order in Q. This means that we can 



simply set e'^^''^ = 1 in Eq. (4.20). Furthermore, the first three terms in Eq.(4.20) correspond to a renormalisation of 
the wall's mass; this is very small and we neglect it here. The two remaining terms contain a coupling of the phonons 
to Q. However, AI5 ^ as ^3 ^ while ttA in the same limit. Thus, the term in gives rise to higher 

powers of Lo — cxk^ than A4i and it is sufficient to consider only the coupling to AI4. The effective coupling between 
the phonons and the wall is mediated by the Fourier transform of the strain tensor Ui = ifc3U2(fc3,r) (since ^2 = 0) 
and the effective interaction Lagrangian is 

L^nt=QJ2'^ksU2{k3,T) (4.21) 



with the coupling constant 



Ck, ==4i7r— (fc3A)3cosech(7rfc3A/2) (4.22) 
2co 



Next, since the coupling is only to phonons in the X2 direction, and with the restriction ki = k2 = 0, we can use, as 
the new Lagrangian for the environment 

Leriv = PvU2{k3,T)u2{-k3,T) + C44:klu2{k3,T)u2{-k3,T) (4.23) 

The mapping to the Caldeira-Leggett oscillator bath is clearly m^j = Pv/L^ and — cxk^ = {C44/ pyY^'^lk^]. 
The way to treat this problem is now straightforward. A total time derivative first needs to be introduced, in order 
to change the coupling between the wall and the phonons from Qu2 to Qu2- This coupling to the velocity of the 
environment can then be eliminated by a canonical transformation on the environmental. The resulting problem 
corresponds to a system couplcd-by its position to a new set of oscillators, with the same ruk and LUk, but a new 
coupling constant Cfc, = cuk^jCkJ^. The dissipative action is now in the usual Caldeira-Leggett form, with a spectral 
function 



Thus the superOhmic J{lu) ^ u!^ is recovered. This is the spectral function we shall use in the discussion of the 
dissipative effects of 1-phonon processes. We now turn to 2-phonon processes. 



D. 2-Phonon Contributions to Action 

Just as in the case of multiple magnons coupled to a domain walll'i, the extra phase space available for the 
interaction of pairs of phonons with a wall, means that 2-phonon couplings give a very different effect from 1 -phon on 
couplings on the wall dynamics. To start with, the summation over the magnetoelastic indices in Eq. ( 4.12 ) is 
obviously more complicated. Let us first consider a static wall, so that Mi = —M.2, -A^s = 0, and let us use only the 
components i?iii — R222 — R333 of the second order magnetoelastic tensor. The summation is then 

\M2\'^[RablRdel - RablRde2 - Rab2Rdel + Rab2Rde2][Gad{K 'i'^n)GbeO<-' , it^n') + Gae{^, il^n)GbdO<-' , it^n')] (4-25) 



We use the fact that ct < cl, which allows us to keep only the term Gt(p, it^^n) in Eq. (4.11) to simplify this 
expression to 

\M2\^RniGT(k,ioJn)GT{k',iujn')[klk'^^{kfk[^ + kjk'^) + klk'^k'i{kl + kl) + k'iklkl{k'^ + k'i)] (4.26) 

Secondly, one must be careful in performing the frequency summations. The two summations (over n and n') should 
be decomposed as summations over the sum and the difference of the frequencies, that is 

n— — oo m— — oo n h. ni 

-TT^ E ')f[l + n^(c.,) + n^(..,0] , + ^ + Mu^u) - n^M] ^ " ) (4.27) 
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where is the boson occupation number, 



By this separation, two different contributions can be identified. If we compare Eq. (4.27) with Eq. (3.7), it is 



clear t hat t he wall can be considered to be coupled to two different oscillator baths. The first term on the right side 



of Eq.( 4.27 ) corresponds to the simultaneous emission of two phonons followed by their re-absorption at a later time, 
while the second term corresponds to the scattering of two phonons off the wall. This latter process is not allowed 
for 1-phonon processes, due to energy-momentum conservation, but is perfectly possible once we consider 2-phonon 
processes. Of course, the scattering term requires phonons to be already present and is thus non-existent at T = (we 
note that precisely the same discussion applies to itaagnons, for 2-magnon processes-involving a single bulk magnon- 
and 3-magnon processes, involving 2 bulk magnonsoLl; the only complication there was the necessity to include a single 
wall magnon as well). The emission/absorption term can be mapped to a Caldeira-Leggett environment provided that 
we identify u> = ojk +0Jk'- The low energy behaviour, w — > then comes from the limit Uk, tok' — > 0, but looking back 
to Eq. (|4]2|), we see that this will implies superOhmic dissipation with a very high power in to, due to the density 



of states of the phonons going to zero as ujk ^ Q- This term can be completely ignored since it will always be much 
smaller than the 1-phonon contribution. 

The scattering term however is very important as it represents Ohmic dissipation. The identification with a 
Caldeira-Leggett bath requires uj = ujk — ujk'^ so that the low energy properties of the bath come from LUk ^k' ■ 



Now, nothing special happens to Eq. ( 4.25 ) in this limit, so the effective density of states for 2-phonon scattering 
processes is extremely reduced. The resulting dissipation is Ohmic, and although the numerical coefficient in front 
of it may be small, .the fact that it is Ohmic means that it may be the dominant dissipative process, especially for 
coherent processestJCd. 

The temperature dependence of the friction coefficient is determined by the phonon density of states, and is 
thus different depending on the dimensionality of the phonons. We examine 1-, 2- and 3-dimensional phonons and 
concentrate on the scattering term, since the emission/absorption term results in superOhmic dissipation. 



1. 3-Dimensional, 2-Phonon Processes 



It is now straightforward to obtain the friction coefficient, although the calculations arc somewhat complicated by 
the anisotropy of the problem. Keeping in mind the restrictions ki + k'l — Q and fc2 + ^2 = coming from the form 
factor of the magnetisation, it is useful to perform the change of variables 



(j) — tan ^ fci/fc2 

UJ = CT{k - k') - CT \{kl + p'f'^ - (fcf + p2)l/2 



(4.29) 



CT{k + k') 



CT 



{ki+py/' + {k'^+p') 



2N1/2 



After doing the angular integrals, the most important contribution to the effective action (the Ohmic one), can be 
isolated to give a dissipative action of the form 



1/T 



dT 



l/T 



dr' 



dujfi2ph{QiT) - Q{t'),uj) I 



■sinh(cj/2T) 



(4.30) 



where the cutoff is taken as the Debye frequency, defined by lod — CTkjy with kjy ~ l/oo, the inverse of the 
lattice spacing. The detailed information of the dissipation is contained in the function fi2ph{Q{T) — Q{t'),uj). We 
are interested in the limit w ^ of the friction coefficient. After a further change of variables x ~ e/T and 
J/ = r-i(e2 - p2)i/2 obtain 



v{Q{r)~Q{T'),0)^ 
3 



^111 

2 

PyCrp 



'-1^ dx 



1 



x^ sinh-'(x/2) 



dyy{y 



2 4 

-X y 



-x')Ml{2yT/cT) cos[2yT(g(T) ~ Q{t'})/ct] 



(4.31) 
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Provided that ko{Q{T) — Q{t')) <C 1 it is possible to expand the cosine and keep only the quadratic term in Q. We 
are then in presence of conventional Ohmic Caldeira-Leggett dissipation. At temperature such that T <^ 9o, the 
Debye temperature defined as ks^D — fi^^D, the dissipation is then characterised by the friction coefficient 

,3d _3 ^2/^^B^ 



where we have reinstated h and ks to get the correct units, is the area of the wall and r(a;) is the Gamma-function. 
Alternatively, it can be written in term of the Debye temperature Qd sls 



where oq is the lattice spacing of the material. The numerical factor comes from taking fc|) — GTr^/cg. It shouldn't 
of course be taken too literally but nevertheless gives the right order of magnitude of the coefficient. Due to the 
presence of the temperature to the 6*'' power, we do not expect the 2-phonon processes to be relevant at really low 
temperatures, the question being of course: how low is sufficiently low! We will consider this question below. 

Notice that we did not include any terms coupling the phonons to the velocity of the wall, as we did for 1-phonon 
processes. The reasons are twofold. First, it is clear that the fact that we obtain Ohmic dissipation coming from 
2-phonon processes is not related in any way to the form of the coupling between the phonons and the magnetisation. 
Therefore, including the terms in Q would only result in temperature dependent superOhmic dissipation, a process 
much weaker than the 1-phonon terms (remember that an apparent Ohmic dissipation in Q is corresponds effectively 
to superOhmic dissipation). Secondly, the component of the second order magnetoelastic tensor required for such a 
coupling .would be of the form i?555 and this is at least four orders of magnitude smaller than i?iii (at least in bulk 
materialsESH) , and can be ignored. 



2. 2-Dimensional Phonons and Neel Walls 



Since the Neel wall is essentially a 1-dimensional structure, the analysis of 1-phonon dissipation in a thin film will 
be quite similar to the one in the bulk. The only difference is that there is now a large component of the magnetisation 
in the direction of the wall motion. The form factors of the Neel wall are thus 

f(^) - \AB) xAN) _ . .(B) 



A^W^A^f) (4.34) 



Again the contribution to the effective action contains longitudinal terms, and transverse terms coming from the 
velocity of the wall. Due to the structure of a Neel wall, the longitudinal term is now proportional to (A33 — ^31)^, 
which is non-zero even in the case of a cubic structure. The transverse kernel, however, still contains the Walker 
velocity cq and we can assume that it still gives the dominant contribution to dissipation and can be used alone. The 
analysis of the Neel wall with a velocity- coupling to the phonons is exactly analogous to what was done for the Bloch 



wall, and the spectral function is completely equivalent to Eq.(4.24) 



JT{-) = ^^^ALc^^cosech^ f ^) (4.35) 
16 Ps CqC^ \ 2ct J 

where Aw is the cross-sectional area of the wall, ps is the surface density and the constant 71.44 is in units of J/ m^ . 

The procedure for analysing 2-phonon processes is similar to the 3-dimensional case, but there are two differences. 
First, the density of states of the phonons is reduced and this will reduce the the power of the temperature in the 
friction coefficient. Secondly, the summation over the components of the magnetoelastic tensor is 

\M3\^[RablRdel - RablRde2 - Rab2Rdel + Rab2Rde2][GadiK 'i'^n)GbeO<-' , il^n') + Gaei^, il^n)GbdO<-' , iUJn')] (4.36) 

If we assume i?iii = R222, then, due to the restriction ki -\- k'^ = coming from the (5— function of the magnetisation 
form factor, the whole summation is zero and there will be no dissipation coming from 2-phonon processes. However, 
there is no particular ground for such an assumption, as the second-order constants are almost unknown. In any case, 
if this is true, then one will certainly couple to the other components of the tensor, say R112, R122 and so on. We will 
assume that such a coupling exists, so that we can write the friction coefficient for the motion of a domain wall as 
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«' ■ mn (M)' ^ m,,y'^n (4^) (J-)' 

where both (R) and p^c^ are in units of Jm~^. 

3. 1-Dimensional Phonons in Magnetic Wires 

The 1-dimensional case is also straightforward to treat. Only longitudinal phonons are present so the relevant 
components of the first- and second-order magnetoelastic tensors are A^a and R^sa respectively. Without specifying 
the magnetisation profile, we simply assume that there exists an average interaction coupling {A)X and ~ (-R)A 
between the wall and the phonons. The 1-phonon spectral function is then 

- --Ji(A>^c.^cosech^ ( ^) (4.38) 



where we assumed that the dominant coupling was still to the velocity of the wall. In cases where such a coupling is 
not present, the Walker velocity cq should be replaced by the longitudinal sound velocity c^. The 2-phonon processes 
still give rise to Ohmic dissipation, with a friction coefficient 

with the expected reduction in the temperature dependence due to the reduction in the density of states. 

E. Ohmic Coupling to Wall Chirality 

The model of chirality tunneling, as proposed by Braun and LossS and Takagi and Tataral^ was discussed in 
Section II. The models considered by these two groups are not completely equivalent. Braun and Loss consider a Neel 
wall in a wire while Takagi and Tatara use a Bloch wall. For definiteness, we concentrate on the Bloch wall model. 
The analysis can be trivially extended to a Neel wall. 

In the model, the magnetisation components are defined as (rni,'m2,m^) = (cos 0, sin ^ sin 0, sin cos 0) with mi 
along the easy axis and normal to the plane of the wall. The angle 9 = 6{x3 — Q) describes a Bloch wall with 
centre at a position Q along the xa-axis. However, the wall is pinned and the dynamical variable is 0(t), with a 
tunneling event corresponding to a transition from (j) = ±7r/2 to = ^t:/2. 

We n ow co nsider the effect of phonons on this chirality tunneling. The magnetoelastic interaction keeps the form 
in Eq. ( |3.18 ), but two points must be noticed: 



(i) The position of the wall is time-independent 

(ii) The form factors of the magnetisation are now explicitly time-dependent, through i^(t). The relevant form 
factors are 

Miiq) = -i(2^)3,5(gi)5(q2)g3A'cosech(^(73A/2) sin^ ^{t) (4.40) 

fe,T) = i(2^)3,5(gi)5(q2)g3A2cosech(^(73A/2)cos2(/.(T) (4.41) 

(93,t) = i(2^)35(gi)5(q2)g3A'cosech(^q3A/2)cos0(r)sin0(T) (4.42) 

We ignore Aif , since it does not depend on the angle 4>{t), and although A4f and ar e non-zero, their effect 



is less important, since they go to zero as —>■ 0, as discussed in Section IV-A (cf. Eq. ([4.9D). There are of course 
no terms depending on the velocity of the wall, and the 1-dimensional character of the Bloch wall still restricts the 
momentum of the phonons in the plane of the wall to be zero. 

It is now straightforward to integrate out the phonon field and to obtain the dissipative action for the chirality 
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^-ff = jz dr dr'TY, d^e^-^^-^ ) ^<i>a(T)$,(r) (4.43) 



47r 



"'0 „ "'O 



with 4'o{t) contained in the functions <i>a(T), with the values of a and h restricted to a, 6 = 2,3,4, and $2 — sin^i/io, 
$3 = cos ^ 00 and $4 = cos sin . This action is not a simple functional of the tunneling variable, but it is 
nevertheless clear that the dissipation is Ohmic, with a dimensionless coefficient a-^ = rj^ /fi given by 

a-^{A)-^^ (4.44) 



Pv 



If we consider a wall of surface — SA^ in Ni (A ~ 500A, Aab ^ IQ^J/nv" and cl ^ lO^m/s), then we obtain 
Q/X ^ 105. Thus, even for a strictly l-dimensional wall, in a sample with all nuclear spins removed, the coupling of 
1-phonon processes to the chirality is obviously going to have a very severe effect on tunneling, and even more on 
coherence. We discuss how severe this effect is in the next section. 

F. Remarks on Power-Counting 

It is possible to explain the frequency dependence of the various spectral functions obtained in this section by some 
simple power counting arguments. 

Let us start with 1-phonon processes. The spectral function describing the effect of isotropic acoustic phonons on 
the tunneling of defects or interstitials in an insulator is of the form J{ijj) ^ uP'^'^ where d is the dimensionality of 
the system. This can be understood in the following way. The square of the matrix element coupling the phonons 
to the defect brings a power or w. The coupling is through a term of the form exp(ik • q(i)) where fc is the phonon's 
momentum and q(€) is the position of the defect. Upon integration of the phonon field, this term results in a factor 
cos(|k-q(i)|) (cf Eq. ( [4.17 ) and ( 4.1^ )) and it is clear that the expansion of the cosine to quadratic order in {q{t)~q{t')) 



results in 2 more powers oi uj — c|k| in the spectral function. Finally, the density of states of the phonons brings a 
factor uj'^~^. 

With the magnetoelastic interaction, the matrix elements are the magnetoelastic constants, but since the interaction 
is through the strain tensor of the phonon field, this will also bring a power of uj to the spectral function. Thus, the 
spectral function of the phonons in a displacement tunneling situation is clearly the results of the uni-dimensional 
character of the Bloch wall, which reduces the effective density of states of the phonons from lu'^^^ to oj, and this, 
independently of the physical dimension of the substrate. 

Coming now to chirality tunneling, it is now obvious that the spectral function relating to the chirality tunneling 
should be Ohmic, since the phonons are effectively l-dimensional, and there is no coupling of the form exp(zA: • (t>(t)). 
Apart from numerical factors of 0(1), the different parameters entering the spectral function are understood easily 
from the spectral fu nctio n of the phonons in a displacement tunneling situation simply by removing two powers of 



k = uj/cl from Eq. ( 4.3^ ) (and with the Walker velocity cq replaced by c^). P- , 

The situation is more complex for 2-phonon processes. It was analysed in detail by Kagan and Prokof 'evE2l for the 
tunneling and diffusion of defects in insulators. In d dimensions, it was found that the Ohmic dissipation resulting 
from 2-phonon processes was described by a friction coefficient with temperature dependence rj{T) ^ y2-i-2d (-^^j^j^ 
"transport effects" included, as is appropriate for isotropic phonons). Again however, this is the result for a perfectly 
isotropic defect. A Bloch wall in a 3-dimensiona.Ljnaterial imposes ki + k[ = k2 + k2 = 0; the problem is much more 
like the problem of magnons coupled to a wallH'El (except that there one has the added problem of wall magnons). 
This then reduces the phonon density of states by a power of 2, which then results in the temperature dependence 
''l^ph ^ ■ III 2-dimcnsions, the only restriction is fci -f- fc^^ = 0, and as such only one power of the temperature will 
be missing with respect to the fully isotropic case, ie., '72p^' ~ ■ Finally, in 1 dimension, there is no lowering of the 

density of states due to the shape of the wall and we get 'ni"2ph ^ 

Now that the dissipative effects of the phonons are known, we go on to examine the various dynamical processes of 
the wall. 

V. TUNNELING OF DOMAIN WALLS COUPLED TO PHONONS AND NUCLEAR SPINS 

We now come to the main practical results of the present paper. In sections III and IV we have seen that the main 
factors determining the tunneling dynamics of a domain wall at low T (for kT < A, the magnon gap) will be (a) the 



23 



coupling to any defect or other pinning potential (b) the random hyperfine potential coming from nuclear spins, and 
(c) dissipation coming from both phonons and nuclear spins. In the pr esent section, we will consider the case where 
the dissipative effects of nuclear spins are small (ie., A/ <C 1, cf. Eq. (3.33)). This will be the case for tunneling in 
wires of Ni or Fe (and compounds of these) . 

We are then left with the problem of a wall moving through a potential coming from both a pinning potential and 
the random hyperfine potential (recall the latter is not small even if A/ ^ 1), with dissipation coming from phonons 
only. Although the random hyperfine potential varies slowly compared to the "bounce time" for tunneling processes, 
it still changes very quickly compared to any external field sweeps in an experiment. 

At first glance the randomly- varying Uhyp appears to make a quantitative analysis of tunneling impossible. How- 
ever, the situation is not so bad for Ni- and Fe-based materials, where the hyperfine coupling is small, and the 
concentration x of nuclear spins is small. Even though Uhyp{Q] {Ifc}) may significantly perturb V{Q), it varies suffi- 
ciently slowly in time that during a tunneling event, the total potential is essentially static. Since in any experimental 
run V is not necessarily known, the main effect of the UnypiQ) is then be to renormalise ^ to a new value. A 
secondary effect is to slightly distort the barrier shape from a quadratic/cubicshape in a way which can be discussed 



quantitatively, for an ensemble of walls. This must be done in terms of Eq. (3.31), and corresponding higher moments 



This kind of theory works best if the change AU in the potential caused by the hyperfine potential inside the wall 
satisfies |AC/|/V^ ^ 1. This is well-satisfied for most Fe- and Ni-based magnets. For rare earth magnets, |AL/|/V^ may 
be very large indeed, and then an analysis using the quadratic/cubic approximation will be almost meaningless. 

Even though |At/|/V^ ^ 1 for Ni- and Fe-based systems, there are still rare large negative fluctuations of AU in 
time. If the external field is swept slowly, the wall can tunnel long before the "bare" pinning barrier is small. Thus 
the observed tunneling characteristics will depend on the sweep rate, and, at slow sweep rates, the wall will cross over 
"too early" to tunneling, from classical activation, and at too high a temperature. 

A complete theory of the tunneling rate in this dynamic nuclear potential turns out to be pitremely difficult; the 
solution of the dynamics in a static potential is difficult enough, even if one ignores tunnelingES. In the following we 
give an approximate theory, valid for slow sweep rates, and for short Ti. Luckily these conditions are appropriate to 
the Ni experiments. This is described below; but first we discuss the tunneling in the complete absence of nuclear 
spins. The tunneling in this case is a simple problem of dissipative tunneling, entirely solvable with the results of 
section IV; we also use these results to look at chirality tunneling and "Bloch" tunneling of walls. 



A. Wall Tunneling in Isotopically Purified Magnets 



We first consider the problem of wall displacement tunneling, ie., the tunneling of the wall from some pinning centre, 
including the dissipative influence of the phonons on this tunneling. We then consider a system without any pinning 
centres (although possibly with more microscopic defects or impurities), and analyse the environmental influence on 
chirality tunneling and Bloch coherent tunneling. 



1. Displacement Tunneling and Phonon Dissipation 

When the defect potential is not perturbed by the nuclear spins, it is straightforward to obtain the corrections to 
the tunneling exponent due to dissipation, since dissipation is fairly weak. We find that the total tunneling exponent 
B{T) is given by 



B{T) = Bo 



l + 15a|'^(r) + -/3('^ 

TT 




(5.1) 



where Bq is the tunneling exponent in the absence of dissipation, given by Eq. ( 2.12 ). This expression includes 
the effects of both Ohmic and superOhmic dissipation, and is valid for walls in bulk, films and wires, through the 
dimensionality d = 1, 2, 3 of the phonons. It is valid up to terms of order (T/Tq)^, ]p*fpi^d that one is not too close 
to the "crossover region" between thermally activated and quantum barrier traversalZj'O'ESI. The dissipative processes 
entering in the effective action are shown schematically on Fig. (Q). As stated earlier, we only consider the superOhmic 
dissipation coming from 1-phonon processes and the Ohmic dissipation coming from the scattering of a pair of phonons. 
To obtain the result (5T), we first consider the bounce obtained in the absence of environment Q{t) — Qq sech(r2oT/2). 
Insertion of this solution into the effective action then yields the factor of lbaf\T) + (6 / tt) (i^'^K The temperature 
correetieft&ito the exponent are found in the usual way by an expansion of the bounce around its zero temperature 
formOt^F^. Since the Ohmic dissipation due to 2-phonon processes is fairly weak and already strongly temperature 
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dependent, it is not necessary to calculate the temperature correction brought by it. Ohmic dissipation should be taken 
into account in the non-zero temperature form of the bounce which is used to calculate the temperature corrections 
due to 1-phonon processes but again, since the dissipation is weak, these can be calculated, in a first approximation, 
as if only 1-phonon processes were present. The calculation is standard and gives the factor ~{<o/tt'^)I3['^\t /Tof . As 
usual, the effect of dissipation is smaller at non-zero temperatures. 



Thus all phonon effects depend on the size of the two dimensionless parameters in (5J_). Ohmic dissipation comes 



from 2-phonon processes and the strength of this effect is parametrised by the dimensionless coefficient af\T) — 
r72p)j(T)/2Miul7o, where in the mass of the wall, Eq. ( |2.6| ), Qq is the oscillation frequency of the wall in the 



potential well, Eq. (2.13) and ?72n)i' friction coefhcient, is defined by the form of the spectral function J{lo) = rjuj 



■ ■ 2ph^ - . 

and given by Eq. (4.39|), (4.37) and (4.32) for 1-, 2- and 3-dimensional phonons respectively. Using the expressions 



for Myj and f2o: we find 

with e = 1 — H/ He defined by the external magnetic field and the coercive field, Eq.(p^), and (i?)^'^^ and (C)^'^^ are 
the relevant second order magnetoelastic constant and the elastic constant respectively. Thus (R)^'^^ is in units of 
J/m , and (C)^^^ ~ pd{c^'^^)'^ , with c^^^ — cl and c^''^ = ct if d = 2, 3 and where pd is the density appropriate to th e 



dimensionality of the material. The constant Kd includes the numerical factors of the spectral functions Eq. (4.32) 



(4.37) and (4.39). For all practical purposes, we can take the order of magnitude estimate Kd ~ 10^^^+^. Notice that 



the power of the ratio between A and oq is the same in all dimensions. This comes by using a wall of surface S^j = clq 
in d = 1 and = Liao, with Li the cross-section of the wall in d = 2. 

The usual 1-phonon processes give the expected superOhmic dissipation. We assume that the coupling between the 
phonons and t he vel ocity of the wall gives the dominant contribution to the dissipation, so that the spectral function 



is given by Eq.( 4.24 ). The strength of the dissipation is given by the dimensionless parameter Pt = fi^o/Mw (coming 



from J(w) = f3uj^) and can be expressed as 



/3i 



id) _ I f (A) \ r(A)A7j 
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(5.3) 



where c*^^-' = cl and = ct if d = 2,3. Similarly, (C)^'*' = Pdc'''^\ with pd the density appropriate to the 
dimensionality of the material. 



We can now analyse the general feature of Eq. (5.1). The salient feature is that there is a competition between 
1-phonon and 2-phonon processes. The 1-phonon part of dissipation is stronger at T = and decreases with increasing 
temperature while 2-phonon processes cause an increase in dissipation with the temperature. Of course, everything 
depends on the relative strength of the dimensionless couplings at and Pt- We now show that both (3t and at are 
quite small, with at < Pt in general. 

Note first that at is strongly reduced by the factor (r/8£))^+'', and to a lesser extent by the ratio of the second- 
order magnetoelastic constant to the elastic constant. This reduction is however offset by the numerical factor 10^''+^, 
the ratio of the width of the wall to the lattice spacing and by the ratio between the magnetisation and the coercive 
field. Finally, there is a weak increase in dissipation as the height of the tunneling barrier is reduced. The total 
value will however be extremely small. Let us consider some parameters appropriate for Nickel (ie., A/oq ~ 10^ and 
R/C - 10^1 and temperature such that T/Qd ^ 10"^. The, we obtain af - lQ-^-'^{Mal Hc)^/'^^-^/^ . Thus, even 
for very small coercive fields and very small tunneling barriers, the effect of 2-phonons processes can be neglected, 
and this even in wires. However, in the region of T ~ lii', so that T/Qd ~ 10~^, the same estimate, in d = 1 gives 

aj^"* ~ 10^^(Afo/-ffc)^^^e~^^'* and it is not obvious that such a term could be neglected. Since the experiments of 
Giordano are performed above 1 K, the question of the relevance of the 2-phonon term remains open. 

The biggest contribution will thus come from the superOhmic dissipation caused by 1-phonon processes. This 
dissipation is independent of temperature, and decreases with decreasing coercive field and tunneling barrier height. 
Again, using the parameters of Nickel (with the first-order magnetoelastic constant A ~ lO^J/m^), we obtain Pt ~ 
{Hc/MoY/^e^/"^. The values of Pt should thus range between 0.01 — 0.1 in an experiment. This result also shows that 
the exact determination of He becomes quite relevant experimentally, as the value of Pt depends mostly on the ratio 
between He and Mq (at least in this simple minded analysis for Nickel). In any case, even for a very strong coercive 
field, it is unlikely that Pt will become greater than 1 so that dissipation does not affect strongly the tunneling rate. 
The dependence of the tunneling rate on the temperature is certainly the most important effect arising from 1-phonon 
processes. In Fig. (^), we show the ratio B{T)/Bq for a Bloch wall in a typical Nickel sample, using 3-dimensional 
phonons, which allow us to neglect 2-phonon processes. 
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Finally, in YIG, the magnetoelastic constants are a lot smaller ( AjC 10 ^ and RjC 10 ^ and we do not 
expect any strong dissipative effect due to phonons, even in a temperature dependence of the tunneling rate. 



2. Block tunneling and Chirality Tunneling? 

So far we have considered only the simple process of the tunneling dcpinning from a defect potential. However, 
suppose that there are no such pinning potentials. What will then be the motion of the wall? One may consider 
one case, where, in absence of nuclear spins or dissipation, the wall would move quite freely (as in the case A S> oq). 
Alternatively, one may consider a situation where a very thin wall, with A/ao '--^ 1 — 10, tunnels through a periodic 
lattice potential (a variant on this uses an artificially imposed periodic lattice potential, acting on a light wall having 
A/ao > 1). 

Let us first address the proposal that one might see coherent "Bloch band" motion of the walls in a periodic potential. 
Notice first just how restrictive the requirements are on such a process. One requires (a) almost perfect periodicity in 
the potential acting on the wall, such that deviations are much less than the bandwidth Aq for translational tunneling 
(b) extremely weak dissipative effects from phonons or other bosonic environments, and (c) no possibility of any phase 
decoherence - by far the most serious source of which is environmental spins (topological decoherence). 

The requirement of near perfect periodicity in the total static potential wall potential Vw{Q) seems almost impos- 
sible to satisfy, at least in the foreseeable future. The basic, pvrohlcm is that the bandwidth Aq will be extremely small, 
even for very small and light walls. Thus, Braun and LossOO propose the example of a YIG wire with cross-section 
lOOAA, with a wall having A ~ 400A; for this very light wall they estimate Aq ~ SOmif (and a crossover to tunneling 
Tq 50mk). Thus, to satisfy the criterion of near periodicity in the potential, one requires that extrinsic corrections 
to VwiQ) (coming from non-uniformities in the wire cross-section, dislocations, defects, etc), be much smaller than 
SOmK. However for this example the wall surface energy ~ 10"'' J/m^ 7 x lO^^K/m^ = 0.07K/ A^, so this is a very 
tough requirement - the displacement of a single atom in the wire will cause a local change in VwiQ) considerably 
larger than Aq!! In reality it is difficult at present to manufacture magnetic wires with better than a 5-10 per cent 
non-uniformity in wall cross-section, so that the fluctuations in the potential will be over energy scales some 3 orders 
of magnitude larger than this. 

Thus the potential felt by the wall will realistically be like that shown in Fig. ^, with some very small periodic 
component superposed. Problpns like this were considered in great detail in studies of the quantum diffusion of defects 
in solids (particularly muons).EZI The actual motion of the wall will be one of quantum diffusive motion between sites 
(eg., neighbouring potential wells), in which a thermal bath takes up the energy difference between sites (here via 
inelastic emission of phonon pairs, this being the dominant Ohmic process at low T). 

Suppose that at t = the wall is found in a potential well as shown on Fig. (^ . Then, if SQ is the distance between 
this potential well and a nearby well with a lower minimum, the wall will tunnel at a rate given by quantum diffusion 



theory ai 



with 

^2ph{T) = Tr,2j>h{T){5Q{Of (5-5) 

where rj2ph{T) is the 2-phonon Ohmic coupling calculated in section IV-C, and ^ is the energy difference between 
the bottom of the 2 wells; the tunneling matrix element /S.{5{Q)) is that between the 2 wells (notice that A((5(Q)) 
decreases exponentially as 5Q increases, so that in most cases 5{Q) will be the distance between neighbouring potential 
wells of the periodic potential). Without more information about the nature of the background random potential, 
one cannot say much more here except that the wall will, at low T, move quantum diffusively towards the bottom of 
the nearest potential well of this potential. 

Turning now to the possibility of chirality tunneling, one notes again the requirement that dissipation and decoher- 
ence be absent, as well as a complete energetic symmetry between the 2 chiral states. The requirement of energetic 
symmetry is still going to be extremely difficult to satisfy, and we also now have the very strong Ohmic coupling of 



single phonons to the chirality (Eq. (4.44)). From the example given after Eq. (4.44) we see that the dimensionless 
Ohmic cauding 3> 1 in almost all cases. This means, according to the standard results for Ohmic T-independent 
couplingjijS that (a) coherence will be suppressed completely (indeed the chirality variable will be localized at low 
T by the coupling to the phonons!), and (b) that the incoherent tunneling rate between the 2 chiralities will be 
suppressed by many orders of magnitude- so much so that we believe that for any reasonably-sized wall it can be 
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neglected (it is probably more to the point to ask why thermally activated transitions between chirality states has not 
been observed at higher temperatures). 

This summarizes the results for an isotopically purified system. We see that the physics of wall displacement tun- 
neling in an isotopically purified system is only weakly perturbed by phonons, in a way easily evaluated theoretically; 
thus it would be very interesting to see experiments performed in such ideal conditions. We expect the results to be 
very different in the presence of nuclear spins, as we now see. 



B. Wall Tunneling in Dynamic Hyperfine field 

We first analyse the tunneling for a Bloch wall; the results for a Neel wall are then obvious. The quasi-static 
potential experienced by the domain wall will be a sum of defect and random hyperfine fields: 

Vw{Q) = ViQ) + UhypiQ;{Ik}) (5.6) 



where V{Q) is given by Eq. (2.7) and UhypiQ', {Ik}) by Eq. ( 3.30| ). Over time scales ^q^, we see UhypiQ; {Ik}) is 
indeed static. 



In addition, the wall is coupled dynamically to phonons, via Eq. ( 3.16 ), and to the nuclear spins fluctuations, via 



At flrst glance it might seem that in systems like Ni, where the nuclear spins concentration x ~ 0.01, and too = 
28.35 MHz is very small, we might begin by ignoring the hyperfine field in favour of the "bare potential" V{Q). 

1/2 

Consider the ensemble averaged value of AC/12 = ^UiU2 between points on either side of the tunneling barrier, here 
Qi is the entry point for tunneling and Q2 is the exit point; recall (Eq. (2.14)) that Qi — Q2 = Xe^^^, in terms of 
e = 1 — Hp/He- We may estimate that 



\AU\ uqI ( X 



V l^OlgHc \No 
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(5.7) 



where V is the barrier height at an applied field He (cf., Eq. ( ^.10 )). For a Nickel wire, using the numbers in section 



II (cf., discussion following Eq. ( 2.15| ), one finds that |A[/|/V^ ~ 0.01 only. Since in any experiment V is not known, 



the main effect of this small UuypiQ) will be to simply renormalise to a slightly different value, which is not directly 
observable, and which does not apparently change the physics very much (of course the same is not true for rare earth 
magnets, where we can easily verify the typically |AC/|/y 3> 1). 

However even for Ni wires the above argument ignores the fact that if the applied field is static, or being swept 
only very slowly, then the actual tunneling of the wall can be controlled by the relatively rare but strong downward 
(ie., negative) temporal fluctuations in Uhyp{Q\ {Ik}) occurring in the tunneling barrier region. 

A general treatment of this problem is beyond the scope of this paper. Here we will give an approximate treatment. 
Let us first consider the case of static external field (ie., zero sweeping rate), with some fixed small value of e. In this 
case we shall define 3 characteristic time scales for the problem, viz., 

(i) the tunneling traversal or bounce time , which will be very short (typically, ilo > 10^ Hz). 

(ii) The characteristic time '^Tjjif f{e) for a fluctuation in the hyperfine field Uhyp to diffuse into or out of the 
barrier; since the barrier width is ~ Ae^^^, we have 

^Tmff{e) ^ ^2^6 > T2 (5.8) 

(iii) Finally, the experimental tunneling time r, ie., the inverse of the tunneling rate F at a given e, in the presence 
of the fluctuating hyperflne potential; we shall estimate this below. 

Let us consider the regime where 

f^o(e)»T2 (5.9) 



^TD^ff{e)-»T-\e) (5.10) 



^-l 



The first inequality simply demands that the hyperfine potential be static during tunneling traversal - since T2 • - 
10'' — 10^ Hz, this is always satisfied. The second inequality corresponds to the assumption that the wall always has 
enough time to sample enough tunneling paths (made available to it by fluctuations in Uhyp through the space of 



27 



possible potentials), that it will sample nearly "optimal" realisations of Uhyp (ie. those giving the maximum tunneling 
r ate , or ne ar to it). We shall see below when this assumption is consistent with our final answer for T"^(e). Once Eq. 
(5.9) and ( ^.10 ) are obeyed, we can estimate T~^(e) by functionally averaging over possible realisations of Uhyp, with 
a tunneling action appropriate to each, to get 



(e) ~ n„{e) / I?[/(g)7'({C/(Q)})e-^»^"'-(«»^^) 



(5.11) 



where 'P{{U{Q)}) is the probability for a particular realisation U{Q) of the hyperfine potential to occur, and 
S(\Vw{Q)}, Q) is the action for the tunnehng trajectory through the total potential Vw{Q) ~ V {Q) + U{Q) (cf., Eq. 
(3.30)). This equation makes no attempts to get pre-exponential factors right. Notice that Eq. ( 5.11| ) is independent 
of T2, because of condition ( 5.10 ). Thus the essential simplifying assumption in Eq. ( ^.ll| ) is simply that if T(e) is 
sufficiently long compared to -^Tjjif f{e), we do not have to worry about the "waiting time" required for U{Q) to first 
approach an optimal reali sati on, every thing boils down to the measure of such optimal configurations. 

The assumptions Eq. (5^) - ( 5.1C ) are probably satisfied in most wall tunneling problems. However, even with 
these assumptions, Eq. (5.11) is not easy to evaluate, it is a complicate d pat h integral over different "paths" U{Q) 
between Qi and Q2- We intend elsewhere to give a thorough analysis of ( 5.11 ). Here we shall estimate the change in 
T(e) brought about by the fluctuating Uhyp{Q) in the weak-coupling regime, where \AU\/V{e) <C 1, appropriate to 
Ni- and Fe-based systems, and show that the change is not important enough to explain the crossover temperature 
discrepancy reported by Giordano et al. 



To make this estimate, we introduce a "typical" potential Ua{Q), having the form 



(5.12) 



for < Q < Qo, where the bare pinning potential V{Q) has a minimum at Q = Qoi and the tunneling end-point is 
at Q = Qo- The value a = 1 refers to the ensemble-averaged "gaussian half- width" value, ie., Eq = |AJ7i2|, so that 



(5.13) 



where N{e) ^ N{X)e^^^ is the number of nuc lear spins swept by the wall in going from Q = to Q = Qq, and 
N{X) — a;(5toA/ao) as before. The ansatz Eq. ( 5.12 ) allows us to reduce ( 5.1l| ) to a 1-dimensional integral: 
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(5.14) 



and where the total potential V^{Q) i s giv en by the sum of V{Q) (evaluated in quadratic/cubic approximation) and 
UhypiQ) (approximated by Ua{Q) in ( 5.12 )), ie.. 



aEolxl^/^ 



(5.15) 



where x — Q /Qo 



Eq. (5.14) is a drastic simplification of (5.11), but it allows us to extract the exponent in the tunneling rate in 
the limit where Eq <C V{e), so that the last term in V^{Q) is a small perturbation. If we write the bare tunneling 
amphtude exponent Bo{e) ~ {MwV{e) Y^'^ (this is equivalent to Eq. ( ^.12 ) or Eq. ( 2.15 ), within a constant ~ 0(1)), 
then we can write the exponential in (5.14) as exp[—B(e,a)/h], where 



B(e,a) =Bo(e) 1 + U^- 



Eq 



V{e) 



2/3 



^3 Eo 

siqn a H — a— 

2 V{e) 



0{{aEo/Vr) 



(5.16) 



where a; > 1. For Ea/V{e) <^ 1, the second term in Eq. ( ^.16 ) dominates over the third one in the integration in 
( 5.14 ) . Doing this integral by steepest descents yields a new tunneling rate 



ro(e)e 



(5.17) 



B(e) =Bo(e)-AB(e) 

^Bo{e)-c{Eo/V{e))B'J'ie) 



(5.18) 
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where ro(e) is the bare tunnehng rate (Eq. ( |2.12| )), and c = ^3/2- 2/9 0.6. Notice that whereas Bo{e) e^/"* (Eq. 
( ^.15 )), we have Ai3(e) ^ e^^^. The decrease in the relaxation time r(e) = f~^(e) is entirely due to the occasional 
downward fluctuations in Uhyp{Q)- However, we emphasise that in the regime of validity of our estimation of Ai3(e), 
ie., for £'o/V^(e) ^ 1, we expect AB{e)/Bo{e) ^ 1 as well (for values of B{e) corresponding to observable relaxation 
times). In particular, for Ni-based magnetic wires, we do not expect a large effect on T(e). 

This concludes our analysis of the influence of nuclear spins on wall tunneling. We conclude that in Ni- and Fe-based 
systems, their influence is fairly small. However, in rare earth systems, where EQ/V{e) S> 1, it is clear that their 
effects on both tunneling and the classical wall dynamics will be very large, and wil l req uire a more sophisticated 
theory, starting from the effective interaction Hamiltonian derived in section III (Eq. ( 3.28| ) and ( 3.29| ). 

Finally, we note that the nuclear spins will also adversely affect any Bloch tunneling, since Uhyp{Q) further disturbs 
any periodicity in the potential. For Ni- and Fe-based systems, this effect is certainly much smaller that that caused 
by disorder and non-uniformities in the wire cross-section; but for rare earth magnetic wires it could easily be larger. 



VI. DISCUSSION: RELATION TO EXPERIMENTS 



Let us now put together what we have learned from the results in sections III-V. With the assumption that for 
temperatures well below the magnon gap energy, the only significant environmental effects on domain wall tunneling 
can come from phonons and nuclear spins, we have set up a fairly complete theory for these effects in the limit where 
the nuclear spin effects arc weak (the limit appropriate to Ni- and Fe-based magnets). We find that in this limit the 
phonon effects can be evaluated more or less exactly at these low energies, and involve 3 main effects - a superOhmic 
1-phonon coupling to the wall velocity, a strongly T-dependent Ohmic coupling to the wall position (both of which are 
usually quite weak) , and a very strong Ohmic 1-phonon coupling to the wall chirality. For a system isotopically purified 
of nuclear spins, these calculations suffice to determine completely the low-T tunneling behaviour; the "displacement 
tunneling" of the wall is weakly modified at finite T from that originally calculated by StampErB. On the other hand 
chirality tunneling is suppressed, and both it and "Bloch tunneling" seem to be practically unobservable, even for 
extremely small domain walls. Addition of nuclear spins makes the problem much harder, but for Ni- and Fe- based 
sytems the effect of the dynamic fluctuating longitudinal hyperfine field on wall displacement tunneling is small enough 
that the change in the tunneling exponent can be estimated. One finds an increase in the tunneling rate over that 
expected without nuclear spins, caused by ocasional strong negative fluctuations in the hyperflne potential. 

The experimental implications for isotopically purifled systems have already been discussed in section V.A. We 
summarize the key result in Fig. (^ , which shows a typical T-dependence of the rate exponent for wall displacement 
tunneling from a pinning potential. It would be of some interest to test this result experimentally on, eg., an 
isotopically purifled Ni magnetic wire, since the theory contains no adjustable parameters once the small oscillation 
frequency fig, the wall coercive fleld _ffc, and the magnetostrictive contants are known. 

For the case where the system contains naturally-occurring nuclear spins, a large number of experimental results exist 
already. The work by Giordano et al. on nominally single domain walls in Ni wires, includes (a) measurements of the 
tunneling rate F, at low T, and the quantum/classical crossover temperature Tq, on various wiresQ (b) measurements 
of the escape fleld statistics over a range of temperatureaJ, and (c) microwave measurements over a-pange of applied 
flelds and temperatures, which give some evidence for level quantization of the trapped wall statea3. In addition to 
these resviits-tliere are many experiments on more complex wall systems, which give various kinds of evidence for 
tunnelingil^'Plr^. These experiments have been done at a variety of sweep rates, and for a wide variety of materials. 

We may make a number of remarks concerning these results, on the basis of the present theory, as follows: 

(i) Quantum/ Classical Crossover. Itis-noticeable that some of these experiments show a crossover temperature Tq 
considerably higher than the predictedEfLl value. In the experiments Tq has usually been determined by looking for a 
low-T crossover to a plateau, in either the relaxation rate or the escape fleld; more sophisticated analyses look at the 
distribution of escape fields. One typically expects Tg to decrease with increasing dissipative coupling to an oscillator 
bath, although how it does so, and the actual T-dependence of the transition rate, obviously depend on the nature of 
the bath (we are at a loss to explain why some experimentalists insist on using formulas appropriate to Ohmic baths 
to describe the crossover; this is never correct unless the low-T dissipation is dominated by electrons, which probably 
only happens for extremely thin walls, if at all). For isotopically purified systems it is clear from Fig. (^), that the 
simple estimate that Tq ^ (rio/Svr), is reasonably accurate. However adding the nuclear spins makes the crossover 
much more complex. Suppose first the applied field is static (or being swept extremely slowly)- then the estimates 
in the previous section apply. We see immediately that for a given temperature, the distribution P{H) of measured 
escape fields will be now be affected by both quantum/thermal fiuctuations and by the fluctuations of Uhyp{Q) in 
the barrier region, ie., there will be a T-independent extra spread in P{H). A rough estimate of this extra spread is 
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obtained from the curvature of -B(e, a) about its minimurp. ia a; this gives an extra spread ~ Eq in bias e, over and 
above the usual T-dependent quantum/thermal spreadingEJEj. 

This then suggests the enticing possibihty that the anomalously high value found by Giordano et al., for Tq in Ni 
wires, might be explained by the fluctuating Utyp, which allows tunneling at a larger bias e than in the absence of 
nuclear spins. However at the prese nt ti me we do not think this is likely, simply because Eo for this case is too small; 
as noted immediately following eq. ( |5.7| ), Eq is roughly 2 orders of magnitude sjpaller than the typical barrier height. 

(ii) Resonant Absorption: In the same way the microwave absorption resultaS will have a linewidth which must be 
increased by UhypiQ); the fluctuations in Uhyp will cause a fluctuation in the resonant frequency, again of order Eq, 
and again T- independent. It is interesting to note that the linewidths in the Giordano- Hong experiments are rather 
large (much larger than one would expect from purely dissipative broadening), and in fact the linewidth is roughly 
what one would expect from the fluctuations in UhypiQ) foi' these Ni wires. It would be useful to test this explanation 
of the line broadening. One obvious way is to redo the experiments on isotopically purified samples. The effect of the 



fluctuating UhypiQ) should also be discernable as a low-frequency (of order Tmffie); cf eq. (5.8)) noise; ie., the 
resonant line is wandering at a frequency far lower than microwave frequencies (and far lower than the frequencies 
of thermal and quantum fluctuations). Similar problems are encountered in strongly non-linear optical systems, as 
well on experiments in glasses. We emphasize here that a full discussion of the effect of these dynamic hyperfine 
expressions on the experiments is premature; it requires the solution of the dynamics of the wall in a time-fluctuating, 
spatially random potential. As noted above, this problem is highly non-trivial, given the complexity of even the 
static problemO. We regard its solution as an important problem- the fluctuating UhypiQ) is so strong in rare earth 
magnets that it is clear that the domain wall motion will be strongly influenced by it, up to temperatures well above 
the hyperfine energy. At low T, the domain wall will be very strongly pinned by the UhypiQ), so much so that the low 
frequency magnetic noise (Barkhausen noise) in the system will be completely determined by the time fluctuations in 
UhypiQ)- preliminary investigations of this question indicate that at low T, the wall motion will be dominated 
by "Levy flights", with a characteristic anomalous diffusion scaling behaviour predicted for the magnetic noise. 
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FIG. 1. The standard model of a Bloch Wall, with axes labelled as in the text. 
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FIG. 2. The standard model of a Neel Wall, with axes labelled as in the text. 
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FIG. 3. Random potential caused by the nuclear spins when kT S> wo- At any given time the potential has a "random walk" 
form (the variance in spatial fluctuations of the potential increase as the square root of the distance travelled by the wall). The 
potential also slowly fluctuates in time, at a rate determined by nuclear spin diffusion. 
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a) 



b) 




FIG. 4. Processes coming into the dissipative action of the domain wall. 3-a) represent the 1-phonon processes, 3-b) represents 
the simultaneous emission and reabsorbtion of 2 phonons and 3-c) is the scattering process giving rise to Ohmic dissipation 
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FIG. 5. Temperature dependence of the tunneling exponent B{T), in the absence of nuclear spins. We show the result for 
3-dimensional phonons, with 2-phonon processes being irrelevant, and assume a value 6/3t/7r = 0.2. 
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FIG. 6. Distortion of the putative Bloch band, caused by imperfections in the background periodic potential. This picture 
is meaningful if the variation caused by imperfections is over length scales long compared to the lattice spacing. 
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